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OVERVffiW 

Tensile  stress-stram  behavior  has  been  widely  recognized  as  the  most  impoitant 
single  property  for  ceiamic  maoix  composites.  In  this  study,  the  tensile  Mure  maehanign^ 
of  the  ceramic  matrix  composites  has  been  investigated. 

Hrst,  the  dog*bone  tensile  ^>ecimens  made  of  unidirectional  ceramic  matrix 
composite  namely:  Nicalon  Hber  (513040,  p  phase)  and  CAS  n  (CaO-Al203'2Si(>2) 
matrix  with  fiber  volume  fractions  of  30%  and  40%  were  tested  in  room  as  well  as  in 
elevated  temperatures.  All  tests  were  conducted  inside  a  chamber  of  a  scanning  clectcon 
microscope  (SEM).  Therefore,  images  at  various  magnifications  revealing  the  intennediaie 
tensile  damage  events  in  the  composite  ^imen  were  obtained  in-sim  along  with  other 
testing  data.  Since  the  lens  in  the  SEM  is  always  kept  at  least  one  inch  away  from  the 
specimen  during  testing,  with  magnificatmns  ranging  from  50x  to  as  high  as  3000x,  the 
heat  on  the  specimen  won't  damage  the  electron  lens.  This  makes  it  posrible  to  obtain 
images  while  doing  tensile  testing  in  high  temperature.  The  damage  histories  of  the 
specimen  at  different  temperatures  were  first  observed  on  the  monitor  of  the  SEM  and 
then  captured  immediately  dither  by  Polaroid  camera  or  video  copy  processor.  Tensile 
stress-strain  curves  were  constructed  at  different  tempMatarcs  (starling  from  room 
tempcraniie,  then  250  OQ  400  «>C.  600  ®C  and  700  ©C).  All  tests  exhibited  a  non-linear 
stress-strain  behavior.  The  efiect  of  specimen  size  (meaning  specimens  with  different 
thidmess  or  width  in  the  gauge  section)  on  the  failure  behavior  was  it  was 

found  that  whfle  thin  specimens  (about  0.0625  "-0.07  "  in  either  thickness  or  width)  tend 
to  have  a  tail  at  failure,  retaining  some  load  carrying  capahiliiy,  thick  specimens  (greater 
than  0.0725"  on  bodi  tiiickness  and  width)  &il  caiastroi^cally,  breaking  into  two  pieces. 
The  effect  of  temperature  on  tensile  failure  strength  was  found  not  quite  significant  within 
the  temperatures  tested.  Micrographs  taken  at  different  loading  and  tenqieraiuies  revealed 
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the  details  of  damage  process  from  first  matrix  cracking  lo  multiple  matrix  cracking, 
i  sequential  random  fiber  slipping  and  fiber  breaking  to  eventual  composite  failure  by 
sudden  fiber  pull-out  and  ooUecdve  fiber  breaking  somewhere  in  the  gage  section.  Matrix 
crack  density  versus  the  stress  with  different  fiber  volume  fraction  and  at  difieient 
temperatures  was  also  studied.  The  results  indicate  that,  for  specimens  with  30%  fiber 
volume  fraction,  the  matrix  crack  initiation  stress  is  less  than  that  of  40%  fiber  volume 
I  fiiacdon.  And  the  matrix  crack  density  (defined  as  the  number  of  matrix  aacks  per  mm) 

I  for  spedmens  of  40%  fiber  volume  fraction  is  higher  than  that  for  those  with  30%  fiber 

i 

volume  fraction.  It  is  also  noted  that  matrix  crack  initiation  stress  was  lower  in  high 
temperature  than  in  room  temperature  for  specimens  of  both  30%  and  40%  fiber  volume 
I  fraction. 

I 

I  Then,  analytical  iix>dels  based  on  finite  dement  method  and  the  singular  integral 

equation  technique  were  used  to  explain  the  tensile  damage  behavior  of  the  oemmic  tnanix 
I  composites.  The  H-shaped  crack  geometry  was  used  in  the  finite  demoit  model  to  study 
the  inteifiice  progression  and  its  effect  on  the  tensile  damage  behavior.  The  results  from 
finite  dement  model  compare  well  with  the  experimental  results.  The  singular  behavior  at 
the  transverse  and  interface  crack  tips  were  studied  using  singular  integral  equations.  The 
stress  intensity  factors  and  strain  energy  release  rates  were  calculated  for  various  crack 
geometry  and  were  used  to  explain  the  failure  mechanism  of  the  composites.  The  results 
from  die  singular  integral  equation  technique  predict  that  once  the  transverse  matrix 
cracks  are  formed,  they  will  propagate  to  the  fibezAnatiix  interface.  This  behavior 
I  conforms  to  the  observed  behavior.  Both  models  assume  that  the  composite  cemsists  of 

I  equal  spaced  fiber  strips  in  the  matrix  material  and  the  problem  is  simplified  as  two 

!  dimensionai. 

I 
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1  INTRODUCTION 

The  leseazch  on  ceiamic  matrix  composite  mateiials  has  intensified  in  recent  years 
due  to  some  of  the  appealing  features  of  ceramics;  such  as  great  stabUiiy  and  resistance  to 
oxidation  under  hostile  (high  temperature  or  corrosive)  environments.  In  comparison  to 
their  metallic  and  polymer  counterparts,  ceramics,  being  brittle  and  low  in  tensile  strength 
and  fiacture  toughness,  traditionally  have  had  little  use  in  structural  applications. 
However,  when  reinforced  with  fibers,  ceramic  matrix  composites  exhibit  an  increase  in 
firacmre  toughness  and  tensile  strength  in  room  as  well  as  high  temperatures  [1-2].  A 
variety  of  ceiamic  matrix  composite  systems  have  been  or  are  being  developed  for 
engineering  applications  ranging  &om  cutting  tools  to  aerospace  structures  [3'4].  For 
examples,  one  can  find  studies  on  the  following  systems  of  ceramic  matrix  composites: 
CVglass  [5],  C/SiC  [6],  SiC/glass  [7],  SiCTLAS  glass  ceramic  [8-13],  SiC/BMAS  glass- 
ceramic  [13],  Siaalumina  [14],  SiC/mullite  [15],  SiC/SiC  (16]. 

In  the  aforemennoned  studies,  experimental  results  on  damage  bdiaviois  of  ceramic 
matrix  con^sitcs  were  obtained  by  failing  the  specimens  with  either  tensile  or  three-point 
bending  loading  at  room  or  elevated  temperature  [9-12].  During  the  theimomechanical 
testing,  the  load  vj.  the  displacement  (and  hence  the  stress  vs.  the  strain)  curves  were 
recorded  but  only  the  postmortem  damage  patterns  were  identified  by  either  a  scanning 
electron  or  an  optical  microscope.  It  is  well  known  that  the  stress- strain  reladons  of  most 
ceiamic  matrix  composites  under  thetmomechanical  loading  usually  exhibit  nonlinear 
behaviors.  Fot  such  ceranric-matrix  composites,  the  failure  mechanisms  that  cause  this 
nonlinear  stress- strain  relationship  is  m^e  complex  than  that  of  their  monolidiic 
counterparts.  They  come  fiom  the  results  of  multiple  matrix  cracking  and  sequential  fiber 
breaking  due  to  weak  inteifacial  bonding  between  the  fibers  and  the  matrix.  The 
approaches  described  above  although  proved  important  data  to  the  overall  understanding 
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of  the  faUiire  mechanisms  of  the  ceramic  matrix  composites,  they  fail  to  correlate  the 
nonlinear  stress-strain  (and  hence  the  sdffoess  reduction)  behavior  with  the  inreimcdiaic 
damage  progression  events  inside  the  cciamic  matrix  composite  specimen.  A  good 
experimental  approach  should  provide  uifonnation  on  the  damage  history  of  the  specimen 
that  can  be  used  to  correlate  the  nonlinear  stress-strain  (and  hence  the  stiffness  reduction) 
behavior  witfi  the  intennediate  causative  damage  events  that  occurred  on  the  ceramic 
matrix  composite  specimen.  Recently,  in  [17-18],  the  tensile  behavior  of  a  hTicalon/CAS  H 
system  was  studied.  Damage  panems  were  identified  and  micrographs  were  taken  to 
capture  the  matrix  crack  propagation.  However,  these  micrographs  were  not  taken  at  the 
same  location  and  therefore  cannot  truly  correlate  the  damage  progression  with  the 
nonlinear  stress- strain  behavior.  And  those  snidies  dealt  only  with  room  temperature.  In 
this  study,  experiments  were  conducted  inside  the  chamber  of  a  scanning  electron 
microsc^e  equipped  with  a  custom  designed  tensile/heating  substage.  This  mffdf  it 
possible  to  directly  observe  and  record  in  situ  the  progressive  tensile  damage  behavior  of 
the  ceramic  matrix  composites  from  the  very  first  matrix  crack  to  complete  fracture  of  the 
specimen  at  any  location  in  the  gage  section.  Test  results  for  both  room  and  elevated 
temperatures  were  obtained.  One  particular  advantage  of  using  SEM  is  for  high 
temperature  testing.  Since  in  SEM,  unlike  in  optical  microscope,  the  electronic  lens  is 
always  kept  at  least  one  inch  away  from  the  heated  specimen,  images  can  be  obtained  at 
magnifications  as  high  as  3000x.  without  damaging  the  electronic  lens  during  high 
temperature  testing.  This  makes  the  technique  very  appealing  for  high  temperature  testing. 

Among  the  factors  that  affect  the  overall  strength  and  toughness  of  the 
matrix  composites,  the  following  are  generally  considered  the  most  important  (a)  the 
thennoelastic  properties  of  the  constituents  0-c.,  the  matrix  and  the  fibers);  (b)  the  relative 
strength  of  the  interface  between  fibers  and  matrix;  (c)  the  volume  fraction  and 
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arrangement  of  the  fibers;  and  (d)  the  ambient  temperature.  Based  on  the  above 
considerations  and  observations  during  the  test,  two  analytical  models  have  been 
developed.  One  model  utilizes  the  finite  element  technique  and  is  based  on  a  periodic 
multiple-row  H-shaped  crack  configuration.  The  other  uses  the  singtilar  integral  equation 
method  with  periodic  transverse  and  interface  crack  configuration.  In  both  models  the  real 
problem  is  formulated  in  two  dimensional  domain. 

2  THE  EXPERIMENTAL  WORK 

The  experimental  work  is  vital  in  this  study.  It  provides  important  insight  and  data 
for  the  understanding  of  the  failure  process  of  this  type  of  composite  material.  And  it  also 
provides  the  foundation  upon  which  analytical  models  are  based.  The  technique  used  in 
this  study  is  believed  to  be  first  developed  in  the  area  of  experimental  mechanics  research. 

2.1  THE  TESTING  PROCEDURE 

Figure  2.1  shows  a  schematic  drawing  of  the  experimental  set-up  which  uses  a 
scanning  electron  microscope  equipped  with  a  tensile  and  heating  substage  to  perform  the 
micromechanical  tensile  testing  of  a  ceramic  matrix  composite  specimen  under  high 
temperature,  fit  this  study,  a  Hitachi  S-2400  scanning  electron  microscope  is  used  which 
is  equipped  with  a  custom  designed  E.F.  Fullam  tensile/heating  combined  substage.  The 
specimen  used  is  simple  “dog-bone”  shaped  specimen  such  as  the  one  shown  in  Fig.  2.2. 
The  actual  experimental  set-up  and  tensile/heating  substage  are  shown  in  Fig.  2.3.  The 
material  used  in  this  study  (Nicalon-fiber/CAS  n  matrix  composite)  was  obtained  from 
Coming  Glass  Works.  Table  2.1  shows  the  thermomechanical  properties  of  the 
constituents  of  the  composite.  Specimens  of  both  30%  and  40%  fiber  volume  fractions 
were  tested.  To  study  the  size  effect,  we  varied  either  the  thickness  or  the  width  of  the 


specimen  for  specimens  with  40%  fiber  volume  fraction.  Fibers  were  unidirectionally 
aligned  with  the  gage-length  direction.  The  tensile  specimens  were  made  by  first  cutting 
the  large  rectangular  ceramic  matrix  composite  panel  as  supplied  by  Coming  Glass  Works 
into  smallCT  rectangular  plates  using  a  Leco  VC-50  Cari/Cut  fine-mesh  diamond  saw. 
Then  the  gage  section  of  the  specimens  was  shaped  using  the  same  diamond  saw  with 
special  a  holder.  The  gage  section  of  the  specimens  was  grounded  to  its  dimensions  using 
a  Dremel  motorized  hand-held  grinder  with  silicon  carbide  and  alumina  oxide  stones. 


Table  1.  Properties  of  Nicalon  fiber  and  CAS  II  matrix 


CAS  n  Matrix 

hKcalon  Fiber 

Nominal  Composition 

Si^C^O 

Elastic  Modulus  Msi(GPa) 

250c 

1 

28.3  (195) 

lOOOOC 

12000c 

11.7  (81) 

22.5  (155) 

Thermal  Expansion(10®/*C) 

5.0  (25®-1000°C) 

4.0  (25“C-1000®C) 

Fracture  Toughness  KjrCMPa  m^^) 

25°C 

2.1610.11 

lOOO^C 

1.3010.17 

Fiber/Matrix  Interfacial  Shear  Strength 

25°C 


15.7±2.0MPa 


7 


* 


To  conduct  high  temperature  testing,  the  bottom  surface  of  the  central  part  of  the 
specimen  as  shown  in  Rg.  2.1,  was  placed  in  direct  contact  with  the  ER  Fullam  heating 
element  (called  heater)  which  is  a  rectangular  plate  made  of  ceramic  material  with  fuse 
wire  circuit  inside.  The  maximum  heated  area  of  the  heating  element  is  0.65"x0.25''.  It 
can  sustain  a  maximum  working  temperature  of  11(X)®C  and  is  equipped  with  a  water- 
cooled  heat  ank  for  continuous  operation.  Temperature  is  measured  by  three  platinum 
30%  rhodium-platinum  6%  iliodium  thermocouples  and  controlled  by  a  stand-along  DC 
power  supply  with  adjustable  voltage  and  current  knobs.  Since  the  whole  operation  is 
conducted  inside  the  chamber  of  a  scanning  electron  microscope  (in  this  study  a  Hitachi  S- 
2400)  which  is  usually  vacuumed  at  1.5x10'^  Pa  or  better,  no  heat-loss  will  occur  due  to 
thermal  convection.  The  top  surface  of  the  specimen  was  first  polished  using  fine-grade 
diamond  paste  containing  15-,  6-,  and  1-micron  particles  until  the  surface  was  well 
finished  and  the  fibers  and  matrix  could  be  seen  clearly  under  a  Nikon  UM-2  microscope. 
Then  the  specimen  was  cleaned  using  a  Bronson  ultrasonic  cleaner  for  10  minutes,  finally 
the  bottom  surface  of  the  specimen  was  coated  with  silver  paint  to  prevent  electric 
charging  and  to  achieve  better  image  before  it  was  mounted  into  the  E.F.  Fullam  tensile 
substage.  Because  it  is  very  difficult  to  drill  holes  in  a  ceramic  composite  spedmen,  the 
top  and  bottom  ends  of  the  specimen  were  mounted  through  stainless  steel  claiiq)S  with 
serrated  teeth  to  the  crossheads  of  the  EJF.  Fullam  tensile  substage.  To  help  in  alignment 
and  prevent  slippage  during  testing  between  the  ceramic  composite  specimen  and  the 
serrated  clamps,  cyanoacrylate-based  extra-strength  epoxy  was  also  applied  on  the  clamp- 
specimen  interfaces.  The  clamped  specimen  was  then  mounted  to  the  tensile  substage  and 
was  left  to  cure  for  at  least  20  hours  to  ensure  that  the  epoxy  had  hardened  and  was 
completely  dry  before  testing. 
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Finally,  the  specimen  together  with  the  tensile/headng  assembly  was  placed  inside 
the  chamber  of  a  Hitachi  S-2400  scanning  electron  microscope  which  is  also  equipped 
with  a  backscatter  detector  to  enhance  the  image.  By  using  the  X-Y  staging  control  of  the 
scanning  electron  microscope,  micrographical  patterns  within  the  central  gage  area  can  be 
observed  and  recorded.  The  tensile  substage  is  driven  by  a  variable-speed  motor  which 
has  a  maximum  speed  of  90  rpm.  Through  a  gear  mechanism  of  100: 1  reduction  ratio,  the 
crosshead  speed  can  be  controlled  within  0.127  mmAnin  (or  0.005  in/min).  The  tensile 
stage  can  provide  a  tensile  load  of  up  to  455  kgs  (1000  lbs).  The  applied  load  was 
increased  gradually  until  the  specimen  failed  totally.  The  damage  to  the  composite  as  the 
load  increased  was  first  observed  on  the  monitor  of  the  SEM.  Then  a  sequence  of 
micrographs  were  taken  to  get  hard  copies  of  various  microcracking  and  damage  patterns 
of  the  specimen  at  different  loading  levels.  Since  the  specimen  is  very  thin  thickness 
(ranges  fix)m  1.0  mm  to  2.0  mm  or  0.04"  to  .08"),  it  is  expected  that  failure  will  occur 
through  the  thickness.  Thus  the  recorded  micrographs  of  the  surface  firacture  can 
represent  through-the-thickness  failure  of  the  specimen.  Quality  of  the  image  and  hence 
the  quality  of  the  micrographs  was  fimher  enhanced  by  transmitting  the  image  signal  to  a 
computer  which  is  equipped  with  an  Imaging  Technology  Advanced  Rame  Grabber 
(AFG)  digital  image  analyze  hardware.  WTtii  the  help  of  the  installed  software,  sharper 
images  were  obtained  through  the  contrast  and  edge  enhancement  operations.  And  text 
can  be  added  to  the  micrographs. 

The  tensile  stage  is  designed  in  such  a  way  that  when  the  load  increases,  the  top  and 
bottom  crossheads  move  in  opposite  directions  to  minimize  the  shift  of  the  observed  site. 
This  is  achieved  by  machining  the  stainless-steel  loading  columns  into  worms  of  reverse 
directions.  Thus  searching  and  refocusing  the  damaged  zone  after  each  load  increment  are 
very  handy.  The  applied  loads  were  recorded  by  a  miniature  load  cell  equipped  with  a 
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digital  readout.  The  relative  displacements  of  the  crossheads  were  measured  by  a  high- 
precision,  strain-gage-type  extcnsometer  which  is  also  fitted  with  a  digital  readout  The 
load-displacement  data  were  recorded  and  converted  into  a  stress-strain  curve.  The 
evolution  of  damage  recorded  by  micrographs  were  identified  with  the  corresponding 
stress  and  strain.  Fmally  the  ruptured  specimens  were  observed  under  a  Nikon  UM-2 
universal  measuring  microscope  for  further  postmortem  examination. 

22  EXPERIMENTAL  RESULTS 
2.2.1  Stress-Strain  Relations 

Typical  stress-strain  curves  for  the  room-temperature,  tensile  damage  behavior  of  a 
NicalonAIIAS  n  with  30%  and  40%  fiber  volume  fractions  are  shown  in  Hg.  2.4  and  Hg 
2.5  respectively.  As  depicted  in  these  two  figures,  the  tensile  damage  behavior  of  the 
ceramic  matrix  composite  specimen  is  characterized  by  a  nonlinear  curve  made  up  of 
three  sections.  The  characteristics  are  the  same  for  all  tests  conducted  in  this  study,  even 
though  microscopically  these  specimens  might  look  quite  different  Upon  load  application, 
the  relation  between  stress  and  strain  was  linear  and  its  slope  was  equivalent  to  the 
stiffness  of  an  intact  Nicalon/CAS  H  specimen  (18.15  Msi  or  125  GPa  for  Vf  =  30%  and 
19.6  Msi  or  135  GPa  for  Vf  =  40%  ).  For  specimens  with  30%  fiber  volume  fraction,  the 
slope  changed  at  about  25  ksi  (Point  A  in  Hg  2.4).  While  for  40%  fiber  volume  fraction, 
the  slope  changed  at  about  30  ksi  (Point  A  in  Hg.  2.5).  First  matrix  crack  is  beUeved  to 
start  at  or  slightly  below  of  point  A  in  both  Vf  =  30%  and  40%  cases.  This  can  be  further 
inferred  from  the  relationship  between  matrix  crack  density  and  tensile  stress  as  will  be 
discussed  later  in  this  section.  In  most  cases,  matrix  cracking  initiates  either  at  the  edge  or 
at  a  location  where  the  fiber  spacing  is  maximum.  Fig  2.6  shows  some  micrographs  taken 
during  one  test  The  second  micrograph  in  Fig.  2.6  clearly  indicates  that  matrix  crack 
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initiated  not  from  the  voids  but  fix)m  the  location  where  spacing  between  the  two  adjacent 
fibers  is  the  largest.  While  in  Hg.  2,7,  matrix  cracks  initiated  fiom  both  the  edge  and  a 
location  where  fiber  spacing  is  very  large  if  not  the  largest  The  matrix  cracks  then 
propagate  perpendicular  to  the  fiber  direction  throughout  the  whole  width  of  the  gage 
section  to  form  multiple  matrix  cracks.  Unlike  the  monolithic  materials  where  failure  is 
controlled  by  a  critical  crack  size,  the  failure  of  the  composite  goes  through  a  process  of 
damage  accumulation.  The  composite  is  insensitive  to  the  voids  and  can  tolerate  very  large 
cracks  before  failure.  At  point  B  (35  ksi  for  Vf  =  30%  and  40  ksi  for  Vf  =  40%),  the 
development  of  multiple  matrix  cracks  reached  a  saturated  stage.  Regularly  spaced  matrix 
cracks  were  formed  in  the  whole  gage  section  of  the  specimen.  This  happened  with  only  a 
small  increment  of  tensile  stress  (about  10  ksi,  as  can  be  seen  fixjm  point  A  to  point  B  in 
both  Fig.  2.4  and  Hg.  2.5).  Further  increase  of  tensile  stress  creates  no  or  very  little 
additional  matrix  cracking.  Matrix  crack  opening,  fiber  debonding,  breaking  and  slipping 
will  dominate  the  rest  of  the  failure  process.  Upon  reaching  point  C  ( about  62  ksi  for  Vf = 
30%  and  58  ksi  for  Vf  =  40%  ),  one  surface  of  the  matrix  cracks  in  the  gage  area  started 
to  open  up  with  crack  opening  displacement  far  more  larger  than  the  rest  of  the  matrix 
cracks.  This  then  was  accompanied  by  massive  fiber  breaking  and  pull-out  in  that  surface. 
And  the  eventual  separation  of  the  specimen  caused  the  load  to  drop  substantially. 
Depending  on  whether  the  specimen  is  thin  (either  in  thickness  or  in  width)  or  thick,  there 
might  be  a  tail  in  the  stress-strain  curve  at  the  load  drop.  Thin  specimens  tend  to  have  a 
tail  at  failure  as  indicated  in  Hg.  2.5  where  both  thin  and  thick  specimens  were  tested.  One 
explanation  is  that  thin  specimens  are  prone  to  bending  during  the  test  This  might  be 
responsible  for  the  lower  failure  strength  and  a  tail  at  failure  for  thin  specimens.  For  thick 
specimens,  failure  is  always  catastrophic.  Fig.  2.8  shows  the  progression  of  damage  for  a 
thin  specimen.  At  failure,  the  specimen  (this  occurred  only  for  thin  specimens)  was  kept  in 
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one  piece  by  the  fibers.  Further  loading  after  failure  was  possible  as  shown  in  Rg.  2.9. 
Micrographs  in  Hg.  2.10  however,  show  the  catastrophic  failure  of  a  thick  specimen. 
Upon  reaching  the  maximum  stress  value,  the  specimen  fails  catastrophically  with  a  very 
big  crack  opening  at  the  failure  surface.  ObsCTvation  on  the  failure  surfaces  of  the 
specimens  vvith  30%  and  40%  fiber  volume  fiiactions,  indicated  that  there  are  agnificant 
differences  in  both  the  amount  and  the  length  of  fiber  pull-out  between  the  two  specimens. 
For  specimens  with  Vf  =  30%,  both  the  amount  and  the  length  of  fiber  pull-out  at  failure 
surface  are  less  than  those  with  Vf  =  40%  as  can  be  seen  from  Fig.  2.11  and  Fig.  2.12.  The 
lesser  amount  of  fiber  pull-out  and  the  shorter  fiber  pull-out  length  is  an  indication  of 
stronger  interfacial  bonding  strength.  This  explained  why  the  tensile  failure  strength  of  Vf 
=  30%  is  slightly  higher  than  that  of  Vf  =  40%.  The  smoothness  of  fiber  pull-out  surface 
as  shown  in  Fig.  2.13  is  evidence  of  non-chemical  bonding  between  fibers  and  matrix 
which  generally  implies  weak  interfacial  bonding  strength.  The  fiber  pull-out  length  and  its 
amount  also  deviated  substantially  among  the  same  batch  of  specimens  as  shown  in  Rg. 
2.14,  implying  that  the  fiberAnatrix  interfacial  strength  may  vary  for  the  same  batch  of 
ceramic  matrix  composite  specimens.  This  is  also  reflected  by  the  fluctuation  of  the 
maximum  tensile  stresses  (Point  Q  among  these  tests. 

Hgure  2.15  is  a  set  of  micrographs  showing  typical  damage  patterns  of  the  spedmen 
at  different  magnifications  after  the  stress  reached  point  B  in  the  stress-strain  curve.  At 
close-up,  one  can  clearly  identify  the  matrix  crack  opening,  fiber  breaking,  slipping  and 
fiber  bridging  of  the  matrix  cracks.  Also  from  Rg  2.15,  one  can  observe  the  typical  "H" 
shaped  crack  pattern  formed  by  intersection  of  the  transverse  matrix  cracks  with  interface 
cracks  during  the  failure  process.  This  "H"  shaped  crack  configuration  will  be  used  in  the 
analytical  models  discussed  later.  Besides  the  easily  identified  matrix  crack  patterns,  fiber 
breaking  also  exhibits  some  patterns  as  shown  in  Hg.  2.16  and  Hg.  2.17.  In  Fig.  2.16,  one 
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fiber  breaks  at  a  location  away  from  the  matrix  cracks  first,  then  a  ray  of  breaking  fibers  is 
form  along  a  line  slanted  away  from  the  line  of  matrix  cracks.  Fiber  breaking  occurred 
mostly  after  the  stress  reached  point  B,  i.e.  after  the  multiple  matrix  cracks  have  been 
formed.  Micrograph  in  Fig  2.17  gives  another  pattern  of  fiber  breaking  in  which  fiber 
breaking  is  randomly  distributed  throughout  the  whole  gage  section. 

The  stress-strain  relations  of  the  Nicalon/CAS  n  composite  at  high  temperatures 
will  be  discussed  later. 

2.2.2  Matrix  Crack  Density 

Matrix  crack  density  defined  as  the  number  of  matrix  cracks  per  1  mm  length  in 
fiber  direction,  was  used  to  characterize  the  failure  process  of  the  composite.  With  the 
help  of  the  scanning  electron  microscope  equipped  with  a  backscatter  detector,  images 
reflecting  the  tensile  damage  pattern  of  the  composite  from  the  first  matrix  crack  to  the 
eventual  failure  of  the  composite  were  captured.  Fig  2.18  and  Rg.  2.19  are  some 
micrographs  showing  the  development  of  matrix  cracks  with  increasing  stresses  for  Vf  = 
30%  and  Vf  =  40%  respectively.  Based  on  the  number  of  cracks  counted  in  the  frame 
shown  in  these  micrographs  and  the  associated  stress  level,  the  matrix  crack  density  versus 
tensile  stress  curve  was  constructed  as  shown  in  Hg  2.20  and  Hg  2.21  for  30%  and  40% 
fibCT  volume  firactions  respectively.  From  the  matrix  crack  density  curves,  it  is  seen  that 
the  matrix  crack  initiation  stress  is  25  ksi  for  Vf  =  30%  and  30  ksi  for  Vf  =  40%.  These 
results  happen  to  coincide  with  point  A  on  the  stress-strain  curves  as  shown  in  both  Rg 
2.4  and  Fig  2.5. 


2.2.3  Effect  of  Fiber  Volume  Fractions 
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As  be  seen  from  the  stress-strain  and  matrix  crack  density  curves  for 
specimens  of  both  30%  and  40%  fiber  volume  fractions,  the  fiber  volume  fraction  affects 
the  tensile  behavior  substantially.  Matrix  crack  initiation  stress  is  less  for  specimens  with 
less  fiber  volume  fraction  (compare  points  A  in  Rg  2.4  and  Rg.  2.5).  This  result  is  in 
agreement  with  that  given  in  [19].  But  the  tensile  failure  strength  for  Vf  =  30%  is  higher 
than  that  for  Vf  =  40%  (points  C  in  the  stress-strain  curves).  This  might  be  the  result  of  a 
relatively  stronger  interfacial  bonding  strength  for  specimens  with  Vf  =  30%.  The  matrix 
crack  density  was  found  to  increase  with  increasing  fibCT  volume  flection  of  the  qiecimen. 
This  finding  also  agrees  with  that  reported  in  [19]. 

2.2.4  Temperature  Effects 

The  same  specimens  were  tested  in  the  SEM  at  250  ®C,  400  ®C,  600  ®C  and  700 
OC.  Rg.  2,22  and  Rg.  2.23  show  the  stress-strain  curves  obtained  at  higher  temperatures 
for  fiber  volume  fiactions  30%  and  40%  respectively.  Within  the  temperature  range  tested 
(  room  temperature,  250  ^C,  400  OC,  600  and  700  ©C),  no  significant  changes  in 
ultimate  tensile  strength  have  been  observed  in  the  stress-strain  corves.  Fig.  2.24  and  Rg. 
2.25  show  the  temperature  effects  on  the  ultimate  tensOe  failure  strength  of  the  specimens 
for  both  30%  and  40%  fiber  volume  fractions.  There  was  a  slight  increase  in  tensile  failure 
strength  for  Vf  =  30%  at  T  =  250  OC  .  Micrographs  showing  the  final  failure  surfaces  of 
the  specimens  with  Vf  =  30%  at  different  temperatures  are  compared  in  Fig.  2.26.  It  is 
seen  that  both  the  amount  and  the  length  of  fiber  full-out  are  the  smallest  for  T  =  250  ®C. 
This  might  explain  why  the  tensile  strength  is  maximum  at  T  =  250  ^C.  For  Vf  =  40%,  this 
increase  of  tensile  strength  happened  at  T  =  400  «>€.  The  phenomenon  of  sUght  increase  of 
tensile  strength  might  be  the  result  of  release  of  thermal  residual  stress  and  changing  of 
interfacial  bonding  strength  due  to  temperature.  Matrix  crack  density  did  not  show 
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dgnificant  changes  with  temperatures  for  Vf  =  30%.  Micrographs  in  Rg.  2.27  show  die 
development  of  matrix  cracks  at  T  =  600  Hg.  2.28  shows  the  relation  between  the 
matrix  crack  density  and  temperatures  (ranging  from  room  to  600  ®C)  for  Vf  =  30%.  But 
for  Vf  =  40%,  it  was  found  that  matrix  crack  density  decreased  with  increasing 
temperature  as  shown  in  micrographs  of  Fig.  2.29  taken  at  T  =  400  ®C  Fig.  2.30  plots  the 
relation  of  matrix  crack  density  versus  temperature  for  Vf  =  40%.  Matrix  crack  initialion 
stress  also  decreased  slightly  at  high  temperatures  for  Vf  =  40%.  Micrographs  in  Hg  2.31 
and  Fig.  2.32  show  the  fiber  breaking  and  slipping  process  in  high  temperatures.  Since  this 
composite  is  basically  a  glass  based  material,  the  stress-strain  behavior  is  expected  to  be 
quite  different  when  temperature  reaches  over  800  ®C. 

2.2.5  Specimen  Sizing  Efibcts 

To  study  the  size  effect,  specimens  of  40%  fiber  volume  fraction  with  different 
thickness  and  width  combination  were  used  in  the  tests.  It  was  found  that  specimens  with 
either  relatively  thin  thickness  or  width  (about  0.0625~0.07")  tended  to  have  a  quite 
different  failure  pattern  than  that  with  thicker  one  (greater  than  0.0725"  in  both  diickness 
and  width).  For  thin  specimens,  first,  the  tensile  strength  was  quite  lower  than  that  of  thick 
specimens;  second,  when  the  specimen  failed,  it  had  a  tail  in  the  stress-strain  curve 
indicating  that  it  still  retained  some  load  canying  capability  and  the  specimen  was  kept  in 
one  piece  by  the  fibers.  The  thick  specimens  on  the  other  hand  always  failed 
catastrophically.  A  comparison  of  tensile  strength  between  thin  and  thick  specimens 
against  temperature  is  shown  in  Fig  2.33.  The  discrepancy  between  the  two  cases  might  be 
the  result  of  bending  that  might  have  occurred  during  the  test  as  explained  earlier.  Size 
effect  of  the  specimens  on  the  flexural  strength  of  unidirectional  carbon  epoxy  composites 
has  been  reported  in  the  literature  [20-21].  Since  there  are  some  definite  relations  between 
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flexural  strength  and  tensile  strength,  tensile  strength  will  be  affected  by  the  size  of  the 
specimen  as  well.  From  the  standpoint  of  Weibull  failure  theory,  thick  specimens  tend  to 
contain  large  defects,  therefore  they  should  fail  at  a  lower  tensile  loading.  Since  the 
composite  with  a  relatively  weak  interfacial  bonding  strength  is  insensitive  to  the  defects, 
the  opposite  results  were  obtained  here.  The  specimen  didn't  fail  by  one  defect  or  crack 
when  it  teaches  the  critical  size,  but  through  the  damage  accumulation  all  over  the  gage 
section.  For  thin  hTicalon/CAS  n  specimens  with  40%  fiver  volume  fi^ction,  the  stress- 
strain  curves  at  high  temperatures  were  shown  in  Fig.  2.34. 
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3  THEORETICAL  MODELS 

As  explained  earlier,  during  testing,  it  was  found  that,  after  the  initiation  of 
transverse  matrix  cracks  at  point  A  (in  both  Figs.  2.4  and  2.5),  with  increasing  load, 
regularly  spaced  multiple  transverse  cracks  were  formed.  When  the  tensile  stress  reached 
a  certain  value,  almost  no  new  transverse  cracks  were  generated  until  failure.  During  this 
process  further  transverse  matrix  cracks  opening,  fiber  debonding,  breaking  and  sliding 
at  the  fiber/matrix  interface  is  believed  to  dominate  the  failure  process.  This  typical 
failure  feature  is  best  described  by  the  periodic  H-cracks  configuration  shown  in  Fig. 
2.15.  Two  models,  one  using  singular  integral  equation  technique,  and  the  other  using 
finite  element  method,  are  adopted  to  either  explain  the  failure  mechanism  or  simulate 
the  observed  tensile  behavior  of  ceramic  matrix  composites. 

The  problem  of  concern  is  basically  a  three  dimensional  problem  because  of  the 
distribution  of  fibers  in  the  matrix.  However,  if  one  assumes  that  the  fibers  are  made  of 
composite  strips  of  width  2Hi,  the  problem  can  be  treated  in  two  dimensions. 

3.1  Singular  Integral  Equations  Formulation. 

It  has  been  shown  that  the  singular  integral  equation  technique  is  a  very  powerful 
tool  in  dealing  with  crack  problems.  Fig.  3.1  is  a  sketch  of  the  proposed  mathematical 
model.  It  is  assumed  that  fibers  are  equally  spaced  in  a  ceramic  matrix,  and  that  the 
thermomechanical  and  fracture  properties  of  the  fiber  and  the  matrix  are  known.  The 
model  contains  cracks  perpendicular  to  as  well  as  parallel  to  the  fibers  at  the  fiber/matrix 
interface  and  which  are  assumed  to  be  periodic.  By  choosing  different  geometrical 
parameters,  one  can  generate  various  crack  geometries.  For  example:  If  one  sets  ai=0, 
a2=H2  and  bi=0  in  Hg.  3.1,  then  the  model  reduces  to  a  geometry  with  periodic  H-shaped 
cracks  as  shown  in  Fig.  3.2.  This  configuration  closely  resembles  the  actual  cracked 
geometry  observed  in  testing. 
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3.1.1  Formulation  of  the  problem 

3.1.1.1  Equilibrium  equations 

The  equilibrium  equations  expressed  in  terms  of  displacements  for  orthotropic  materials 


are  as  follow: 

32u  3^u  -  3^  ^ 

Pl3x2'*‘3y2‘^P3axay  =  ® 

3V  f.  3^  p  3^ 

9x2  ^  ^9y2  ^^9x9y 

where  Pj,  Pj  P3  are  material  constants  [23] 

Assume  that  the  solutions  are  of  the  following  forms: 
u(x,y)  =  u<^>(x,y)  +  u®(x,y) 
v(x,y)  =  vt^KXfy)  +  vPKx.y) 


where 


u(l)(x,y)  =;  f  f(a,x)cosayda 
%  Jo 

vO)(x,y)=§f  g(a,x)sinayda 
it  Jo 


(3.1a) 

(3.1b) 


(3.2a) 

(3.2b) 


(3.3a) 


(3.3b) 


u(2)(x,y)=“£  h(a,y)smaxda 

v(2)(x,y)=-f  ((a,y)cosaxda 
Tt  Jo 


(3.4a) 


(3.4b) 


here  f(a,x),  g(a,x).  h(a,y)  and  t(a.y)  are  shape  functions 

Substituting  (3.3ab)  and  (3.4ab)  into  (3.1ab).  one  gets  the  following  characteristic 


equation: 


i4  +  p4r2  +  P52  =  0 


where 


and 


p5=Vp2/Pl 


and  the  roots  are  obtained  as: 


where 


ri  =  -  r3  =  Wj  +  iw2  = 


r2  =  -  r4  =  W3  +iw4  =  V(-P4  ’ 


Then  we  have: 

f(a,x)  =  A(a)erica  +  B(a)e-riC«  +C(a)er20a  +  D(a)e-r20«  (3.6a) 

g(a.x)  =  p7[A(a)eriax  .  B(a)e-riax] 

+  p8[C(a)er2«x-D(a)e-r2«x]  (3.6b) 

h(a,y)  =  E(a)e>’iay/P5  +  F(a)e-riay/p5 

+  G(a)er2ay/p5  +  H(a)e-r2«y/p5  (3.6c) 

((a,y)  =  p9[E(a)eriay/p5  -  F(a)e-riay/P5] 

+  Pio[G(a)  er2ay/P5  -  H(a)e-r2ay/p5]  (3.6d) 

Here  two  types  of  material  will  be  distinguished  according  to  whether  the  roots  of 
the  characteristic  equation  are  real  or  complex. 

Material  type  I:  where  both  rj  and  r2  are  real  numbers(w2=W4=0). 

Material  type  H:  where  rj  and  r2  are  complex  numbers. 

Applying  symmetry  conditions: 

u(x,y)  =  -u(-x,y)  (3.7a) 

v(x,y)  =  -v(x,-y)  (3.7b) 

we  obtain: 

B(a)  =  -A(a)  D(a)  =  -C(a)  (3.8ab) 

F(a)  =  -E(a)  H(a)=-G(a)  (3.8cd) 


Material  type  I  will  be  considered  because  most  materials  fall  into  this  category.. 
Based  on  the  fact  that  both  u  and  v  vanish  when  y  -»  00,  and  that  the  problem  of  concern 
is  symmetrical  about  the  x  axis,  the  displacement  functions  can  be  written  as: 


20 


4  f** 

u(x,y)  =  —  J  [A(a)sinh(wicxx)  +  C(a)sinh(w3ax)]cosayda 

+ “  f "  [E(a)e-lwilay/P5  +  G(a)e-lw3lay/p5]sinaxda 
ic  Jo 

4  f®* 

v(x,y)  =  —  I  [P7A(a)cosh(wiax)  +  PgC(a)cosh(w3ax)]smayda 
7C  Jo 

--  f  [p9sign(wi)E(a)e-'wilay/p5  +  piQsign(w3)G(a)e-lw3l«y/P5]cosaxda 
K  Jo 

3.1.1.2  Stress-strain  relations 

Let  A  = 

then  for  orthotropic  materials,  under  plane  stress  conditions' 

'Cxy(x,y)  =  GxyYxy 

where  Ex,  Ey,  Gxy,  Vxy  and  Vyx  are  the  othortropic  material  constants. 

3.1.1.3  Strain-displacement  relations 


In  terms  of  u  and  v,  the  strains  are  given  as: 
3u  ^  j 

^  ~  ax’  ^y  “  ay  ~  ^ay  ax^ 


(3.9a) 


(3.9b) 


(3.10a) 

(3,10b) 

(3.10c) 


(S.llabc) 

Substituting  eqns  (3.9ab)  into  (3.11)  and  then  into  (3.10),  the  stress  can  be  written 


as: 

7C(l-VxvVvx) 


^xy^yx)g^(x  y)  _  f  [YjE(a)e'l'^i'®y/P5  +  Y2G(a)e"lw3lcxy/P5](xcosaxda 
2cx  Jo 

+  [  [2y3A(a)cosh(wiax)  +  2Y4C(a)cosh(w3ax)](xcosayda  (3.12a) 
Jo 

^^^~2^^y^~ay(x,y)  =  I  [Y5E(a)e‘lwilay/p5  +  Y6G(a)e‘'w3lay/P5]acosaxda 
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+  J  [2Y7A(a)cosh(wiax)  +  2YgC(a)cosh(w3Cxx)]otcosayda  (3.12b) 
~^xy(x»y)  =  J"[YiiE(a)e'''^i'“y^5  +  Yi2G(a)e''^3'ay/p5]asinaxda 

+  J  [2Y9A(a)sinh(wiax)  +  2YioC(a)sinh(w3ax)](xsinayda  (3.12c) 

The  above  e;q)ressions  are  for  plane  stress  problems.  It  can  also  be  used  for  plane 
strain  problems  with  the  following  substitutions: 

'’*y=-||Vxy. 

£yA  “  Ex  and  ExA  ®  Ey, 

3.1.2  Boundary  and  continuity  conditions 


At  the  interface(xi=Hi,  X2=-R^ 


^ixx(Eii»y)^2xx(*^2»y) 

0<y<oo 

(3.13a) 

'^ixy(Hi.y)=^xy("K2»y) 

0<y<‘>® 

(3.13b) 

Ui(Hi,y)=U2(-H2,y) 

0<y<bi  or  b2<y<®® 

(3.14a) 

v,(Hi,y)=V2(-H2,y) 

0<y<b|  or  b2<y<®® 

(3.14b) 

<Jixx(Hi.y)=P3(y) 

bi<y<b2 

(3.15a) 

'Cixy(Hi,y)=P4(y) 

bi<y<b2 

(3.15b) 

At  y  =  0  (cracks  normal  to  the  interface) 

Vi(xi,0)=0  -Hi<Xi<-ai  and  ai<X|<H| 

(3.16a) 

V2(x2,0)=0  -H2<X2<-a2  and  a2<X2<H2 

(3.16b) 

ajyy(Xi,0)=-Pi(Xi) 

(3.17a) 

O2yy(X2,0)=-p2(X2) 

-a2<X2<a2 

(3.17b) 

(3.18a) 

‘C2xy(^2>®)^ 

-H2<X2<H2 

(3.18b) 

At  xj  =  0  or  X2  =  0 

Ui(0,y)=0 

0^y:Se>o 

(3.19a) 
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U2(0,y)=0 

0<y^ 

(3.19b) 

‘^ixy(a»y)^ 

O^y^ 

(3.20a) 

't2xy(a*y)^ 

O^y^ 

(3.20b) 

where  p,(X|)  is  the  normal  traction  on  the  fiber  crack  surfaces  (-ai<x,<a,  and  y=0), 
p2(x2)  is  the  normal  traction  on  matrix  crack  surfaces  (-a2<X2<^^  y=0)»  PsCy)  is  the 

normal  traction  on  crack  surfaces  along  the  interface  (X|=Hi,  bi<y<b2)  and  P4(y)  is  the 
shear  traction  on  crack  surfaces  along  the  interface  (Xj=Hi,  bi<y<b2). 


Define: 

avi(x„0)  aviM 

9X2 

<|)3(y)  =  ^Vi(+Hi,y)  -  V2(-H2,y)].  <|)4(y)  =  ^Ui(+Hi.y)  -  U2(-H2,y)]  (3.21abcd) 


Then,  from  (3.14ab)  and  (3.16ab),  one  can  write: 

(|>ifX|)  =  0  aj<lxjl<Hj  (3.22a 

(jijCxj)  =  0  a2<lx2l<H2  (3.22b 

<])j(y)  =  0  0<y<bi  &  b2<y<®®  (3.22c 

(j)^(y)  =  0  0<y<bi&  &  b2<y<«»  (3.22d 

Using  eqn(3.9a)  and  (3.12c),  eqns(3.19a,b)  and  (3.20a, b)  are  satisfied  identically 

Substituting  eqns(3.18a,b)  into  (3.12c)  we  have: 

E(a)  =  -^(a),  E*(a)  =  -3^*(a) 

Yii  Yii 

Then  eqns(3.9a,b)  and  (3.12a,b,c)  become: 


(3.22a) 

(3.22b) 

(3.22c) 

(3.22d) 


u(x,y)  =  —  f  [A(a)sinh(wi  ax)  +  C(a)sinh(w3ax)]cosayda 
nJo 

+-  f~  E(a)[e-'wi'«y/p5-^e->W3lay/P5]sinaxda 

It  Jo  Yi2 

4 

v(x,y)  =  —  [^A(a)cosh(W|ax)  +  |58C(a)cosh(w3ax)]sinayda 
it  Jo 


(3.9a*) 
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f"E(a)[Mgn(wi)e'lwilay/p5.%,osign(w3)e-lw3lay/p5]cosaxd^  (3.9b*) 

JC  Jo  *12 

^  p  E(a)[Yie->Wilay/P5  -  Y^e-lw3«ay/p5]acosaxda 

Jo  il2 

+  J  [2Y3A(a)cosh(Wjax)  +  2Y4C(a)cosh(w3ax)](xcosayda  (3.12a*) 

?t(l-yyx)- ^  ^  f"  E(a)[Y5e-'Wilay/P5  -  Y^e-IW3lay/p5]cxcosaxda 

2Ey  y  Jo  il2 

+  J”  [2Y7A(a)cosh(wiax)  +  2Y8C(a)cosh(w3ca)]acosayda  (3.12b*) 

2^-^y(x,y)  =  J"  E(a)Yii[e-'wi'«y/p5  -  e-lw3lay/p5]asinaxda 
+  J”  [2Y9A(a)sinh(w,ax)  +  2YioC(a)sinh(wjax)]asmayda  (3.12c*) 


Applying  eqns(3.16a,b),  at  y  =  0,  one  obtains: 

Vi(xi,0)= E(a)[p9sign(wi)-^iosign(w3)]cos<xxida  ai<Xi<Hi  (3.23a) 

V2(X2,0)=  I  E*(a)[p9*sign(wi*)-^^^io*sign(w3*)]cos<xx2da  a2<X2<H2  (3.23b) 

with 

Y,3E(a)cismax,do, 

J"  Y,5»E*(a)asmoxjda 
The  inverse  Fourier  transform  gives: 


<l)i(xi)sin(xxidxi 


(3.24a) 
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Substituting  eqns(3.24a,b)  into  eqn(3.9b)  respectively  and  applying 
eqns(3.17a,b),  one  can  derive: 


f*i  r  [2Y7A(a)cosh(Wiaxi)  +  2YgC(a)cosh(w30tXi)3(xda 

Jo 


Laj  t-Xi 


(3.25a) 


Yu*  P  ^^t  +  f**  [2Y7*A*(a)cosh(wi*(xx2)  +  2Y8*C*(a)cosh(w3*ax2)lada 
J-aj  t-X2  Jo 

where  Yi  and  YL*  are  given  in  Appendix  B. 

Substituting  eqns  (3.24a, b)  into  eqns(3.15a,b),  one  has: 
it(l-VxyVyx)  ^  ^  r  f“‘^sinatdt[Yie-'wilay/p5-Y^e-lw3lay/p5]coscaida 

+lim  f"[2Y3A(a)cosh(Wiaxx)  +  2Y4C(a)cosh(w3axi)]acosayda  bi<y<b2  (3.26a) 
x,-4H.  Jo 

:r;J-p.(y)  =  lim  T  p^sinatdtYnle-'^i'^y/Ps  -  e-<W3lay/P5]sinocx,da 
2Gxy  ».-♦**,  Jo  Jo  113 

+lim  f  [2Y9A(a)sinh(W|axi)  +  2YioC(a)sinh(w3(xxi)]asinayda  bi<y<b2  (3.26b) 
Jo 

Applying  eqns(A.5a,b),  eqns(3.26a,b)  can  be  further  reduced  to: 


Ha 


•]dt 


yiCO,  H|-t _  .JiU - - ,1 

2Yi3  ^Sri2 

+lim  r  [2Y3A(a)cosh(wiaxi)  +  2Y4C(a)cosh(w3axi)]acosayda 


bi<y<b2 


(3.26a*) 


and 
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h. 


jwali 

_A_ 


■]dt 


••l  Ps^ 

+lim  P  [2Y9A(a)sinh(WiOCXi)+  2YioC(oc)siiih(w30CXi)]asinayda 

Jo 

bi<y<b2 


Applying  cqns(3.14a.b),  (3.13a.b)  and  (3.21a,b,c.d).  one  obtains: 


A(a)sinh(wiaHi)  +  C(a)sinh(w3CcHi) 

+  A*(a)sinh(Wi*(xH2)  +  C*(a)sinh(w3*aH2)  =  Ri(a) 
p,A(a)cosh(Wi(xHi)  +  p8C(a)cosh(w3(xH,) 

-  p7*A*(a)cosh(wi*(xH2)  -  p8*C*(a)cosh(w3*(xH2)  =  RjCa) 

[Y3  A(a)cosh(Wi(xHj)  +  Y4C(a)cosh(w3(xHi) 

-XiY3*A*(a)cosh(wi*(xH2)  -  XiY4*C*(a)cosh(w3*aH2)]a=R3(a) 
[Y9A(a)sinh(WiaHi)  +  YioC(a)sinh(w3(xHi) 

+  X2Y9*A*(a)sinh(Wi*otH2)  +  X2Yio*C*(a)sinh(w3*aH2)]a  =  R4(a) 

where  Ri(a),  RjCa).  R3(a)  and  R4(a)  are  given  in  Appendix  A. 

Solving  eqns(3.27abcd)  for  A(a).C(a).  A*(a)  and  C*(a),  we  have: 


(3.26b*) 

(3.27a) 

(3.27b) 

(3.27c) 

(3.27d) 


A(“)  =  cosh(B,aH;? 

c*c«) 

where  f(a),  gi(a),  hi(a),  mi(a)  and  ni(a)  (i=l  ,2,3,4)  are  given  in  Appendix  C. 


(3.28a) 

(3.28b) 

(3.28c) 

(3.28d) 


Substituting  (3.28a,b,c,d)  into  (3.25a,b)  and  (3.26a*,b*): 


26 


-ri-" 

n  I  4-x 

J-aj 

/ 


«i  •'*‘2 


•bj  pbj 

+  j  K,3(Xj,t)(J)3(t)dt+  I  Ki4(Xi,t)(l>4(t)dt 

lb,  \ 


r. 


K  I  n-Xj 


K2i(X2,t)(j>i(t)dt  +-  j  [~  +  xK22(X2,t)]<|)2(t)dt 
J-a,  •'-aj 


J.b2  |.bj 

K23(X2.t)<t)3(t)dt+  K24(X2,t)(l)4(t)dt 

b,  \ 

(^~^*yv^*vx) 


xyr  yx^u  (x ) 

2Y,4*Ey* 


J^K3,(y,0Wt)dt  +  | 


K3i(y.t)(l)i(t)dt+  I  K32(y,t)<|)2(t)dt 
-a,  •'•aj 


J»b2  ^bj 

K33’(y.t)<|)3(t)dt  +  K34’(y.t)<l)4(t)dt 

b,  Jb, 

(l-V^yVyx) 


f  K34’(y.t 

Jb. 


2Ev 


bi<y<b2 


r 


K4i(y,t)<l)i(t)dt  +  I  K42(y,t)(l>2(t)dt 

J-a,  Jti2 


+  [  K43’(y.t)(l)3(t)dt  +  [  K44-(y.t)<I)4(t)dt 

Jb,  alb, 


f  K«'(y.< 

Jh. 


{329&) 


(3.29b) 


(3,29c) 
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=  2g  "P4(y)  bi<y<b2  (3.29d) 

xy 

where 

Kn(xi,t)  [kj(x,.a)e-aCHrt)Mwil  +  k2(xi.a)e-a(Hi-0p5/lw3l  ]da 

Ku(xi.t)  =:~jr  0£3(x,.a)c-«(H2-t)p5*/lwi*l  +  k,(xi.a)e-a^^^^^ 

Ki3(xi,t)  =:j“r  [JiCxi.a)  +  JjCxi.a)  ]cosatda 

Ki4(xi,t)  =  [JsCxi.a)  +  J4(xi.a)  ]sinatda  (3.30abcd) 

K2,(x2,t)  [k5(x2,a)e-«(Hrt)Mwil  +  k6(x2.a)e-a(Hi-t)Mw3l  ]da 

K22(x2,t)  =:“J7  |:k,(x2,a)e-a(H2-t)p5*/lwi*l  +  k8(x2.a^^^  ]da 

K23(x2.t)  +  J6(x2.a)  Icosatda 

K24(Xi,t)  =  [J7(x2.a)  +  J8(x2.a)  Isinatda  (3.3 labcd) 


+  lim  -r  [k9(y,xi.a)e-a(Hrt)Mwil  +  kio(y,Xi,a)e-aCHrt)p5/lw3l]da 

TtJO 

Kj2(y.t)=  to  if  [k„(yjt,.a)e-a(H2-t)Pi*/lw,*l  +  k,,(y,x.,o)e-«(H2-0fe^^^^ 

W  =  if 

K^'(y.t)  =  to  if  (3.32abcd) 

lw,lv  Iwtiv 

,  ■  Yn  I- _ §5 _  _ §5 _ 

+  to  ip  [k|3(yjc,.a)e'“(Hr<)p5^w,l  +  k„(y,x,,o)e'“(H,-t)p5flWjl  jao 

x,^M,  ItJO 


28 


K42(y,t)  =  lim  “f  [k,5(y,Xi,a)e'®('^2*t)Ps*/*Wi*l  +  k,5(y,Xi,a)e'®(^2‘t)Ps*/'w3*l  ]da 

*,->11,  TC-'O 

V(y.t)  =fa  if  + 

= j“ir 

The  derivation  of  the  terms  kj  and  Jj  (i=l  to  16)  are  given  in  Appendix  C. 

Note  that  for  t  =  y  the  integrals  (3.32c,d)  and  (3.33c,d)  are  divergent.  These 
divergent  parts  must  be  studied  and  separated  by  analy2dng  the  asymptotic  behavior  of 
the  integrands. 

After  separating  the  singular  parts  in  (3.32c, d)  and  (3.33c, d),  eqns(3.29c,d)  are 
finally  reduced  to: 


ri  r*2  ^**2 

K3i(y,t)<l)i(t)dt+ J  K32(y,t)<l)2(t)dt +^<l)3(y)+  J  K33(y,t)(|)3(t)dt 


I  K:,4(y.0W0dt 


t+y 


bi<y<b2 


(3.29c*) 


K4i(y,t)(l)i(t)dt+  K42(y.t)(})2(t)dt  + 

J-aj  J-aj  .'b, 


S-y  <+y 


JifeWdt 


+  J  K43(y,t)<|)3(t)dt  -  ■^<t)4(y)  +  J  K44(y,t)<l)4(t)dt 


,b2 


=  2^  ■P4(y)  bi<y<b2  (3.29d*) 

xy 

where  p,,  p2,  P3  and  P4  and  K34(y,t),  K43(y,t)  and  K44(y,t)  are  defined  in 

Appendix  D  The  additional  single-valueness  conditions  for  each  crack  configuration  will 
be  discussed  later. 
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3.1.3  Asymptotic  Analysis  of  the  kernels 

Depending  on  the  aack  geometries,  some  or  all  of  the  kernels  might  be 
unbounded  as  a  oo.  The  unbounded  kernels  will  affect  the  values  of  singularities  at  the 
crack  tips.  Therefore  an  asymptotic  analysis  is  necessary  to  determine  these  smgularities. 
Generally  a  kernel  can  be  expressed  as: 

Ky  =  Kijs  +  Kijf  (ij  =  1.2,3.4)  (3.34) 

where  Kjjs  is  the  unbounded  part  and  Kjjf  is  the  bounded  part 

Let  kj^  (Xi.a)  (i=l-8)  be  the  asymptotic  part  of  IqCxi.a),  when  a  then  we 

have: 

K„s(x,.t)  +  k,-  (xj.a).-®'-*’""'  ]da 

+  rcosh(w3axi)e‘“^^‘"‘^^ 

+  rcosh(wjax,)e-^^«‘-‘^‘''"‘'"”‘^‘'’da 

^  rcosh(w30cx,)e"t^-^^-^'"»^”»'"»"da 

1  (H,-t)P,/lwiH-Hilwil 

1C  -  t)p5  /Iw,  l+Hj  I  Wi  -  (WiXi  f 

1  (H,-t)P,/lw.l+Hilw3l  _ 

1C  ^  [(Hi  -  t)p5  /  Iwi  l+Hi  1 W3  -  (W3X1 

1,  (H,-t)Ps/lw^»+Hilwil 

1C  [(Hi  - 1)|35  /IW3  l+Hi  I  WiO^  -  (WiXif 

.U  (Hi-t)p^/lw3H-Hilw3l  (3  35) 

1C  ^*[(Hi  -t)p5  /IW3  l+Hj  IWjl]^  -  (W3Xi)^ 
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Similarly  one  can  derive: 


If  .  ,v  U  (H,-t)p;/lw:i+H.lw.l 

i2sC  i>)  „^3[(H^_t)p;/lw;i+H,lwJl*-(WiXi)* 

+  1^  (H,-t)p;/lw:i+H.lw3l 

+  I2  (H,-t)p;/lw;i+H.lw,l 

[(Hj  -Ops/Iwjl+Hjlwilf  -(WjXi)* 


1  (H,-t)p;/lw;i+H.lw3l 

n  ^  [(Hj  -Opj/Iwjl+Hjlwjlf  -(W3X1)* 

If  rv  =  (H.-t)p3/iw,l4-H3lw;i 

2isC  2.)  „^»0J[(H,-t)p5/lw,l+H3lw;if-(w;xj)* 

+  ii  (H,-t)p3/lw.l+H3lw;i 

n  '“[(Hi - t)Ps/lw,l+H2lw;if -(w;xj)" 

+  I1  (H.-t)p3/lw3l+H,lw:i 

K  '“[(H,  -Dps/Iwjl+Hjlwllf  -(w*xj)" 

+  I1  (H.-t)P3/lw3l+H3lw;i 

7C  '°^[(Hi  -t)p5/IW3l+H2lwJlf  -(W3X2)* 


^22s(*2»*)  “  ^09 

% 


(H2-t)p;/iw:i+H2iw:i 
[(Hj  -t)p;/iwi*i+H2iw;if  -(w.xj)' 


+  I1  (H.-t)p;/iw:i+H2iw;i 

It  “®[(H2  - t)p;/i wi*i+H2iw;if  -(wixj)^ 


+ii  (H2-t)p;/iw;i+H2iw:i 

jt  '”[(H2-t)p;/iw;i+H2iw;i]"-(w;x2)" 


+lx  (H2-t)p;/iw;i+H2iw;i 

Jt  ""[(H2-t)p;/iw3*i+H2iw;i]^-(w;x2)^ 


The  asymptotic  expressions  for  KjiCy.t),  K32(y,t),  K4i(y,t) 


(3.36) 


(3.37) 


(3.37) 

and  K42(y,t)  are 


obtained  as  follow: 
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1  J  J  ^  ^ 

KsisCy.t)  =  2Y„TC*-'^KIw,ly)^  ^ 


H,.t 


i 

-j-m  t^2  'Ti2  v'^3'y^' t\2 

P^2  '*'^1  ft  2  '*‘(^1"^' 


P5^ 


+  lim  “f  [k9.(y.Xi,a)e-*^'-‘^'*‘'^''  +  k,o.  (y.Xi.a)e 

x,-m,  TCJO 

__l_r<, _ !k! _ ^ — ^1 

+  lim  Ya^si  +Y4\g. 

*t“***i  ^^13^80  ^ 

+  lim  Y3^.ia±Y4^ii  rcos(ay)e'®'^‘"‘’'*’'^’'^^‘‘*‘’'*’'’da 

t-**.  irV—l—  *0 


•o(H|-t)P|/lw,l 


]da 


*i**'*i  ^a^go 

1  ,  _ Iki 


2Yi37?  ^\IWily)!^(H^-t)2 


...Illv _ ikl 


Ps* 


+(H,-t)2 


...  (H,-t)Mw,l  _  (H.-t)p5/lw3l 

Ps  r/TT  *\Q  11.,.  n2  •  ,,2  Pfi 


[(H,-t)Ps/IWilf +y*  ■  ''‘’[(Hj-DPs/IWjlf+y* 


^32s(y>^)  “  P? 


(H,-t)p;/iw;i 


+  P8 


(H,-t)p;/iw;i 


[(Hi-t)p;/iw;i]"+y"  ■  ''*[(H2-t)p;/iw;if+y" 
IWilv  Iwdv 

K-  /„rt=Jli-r - fe - . - Si. 

2y„7cWily)V  (Iwi’y)!. 


i 


(3.38) 

(3.39) 


K42s(y»^)  “  *  Pll 


‘’*'[(Hj-t)p5/lw,lf +  y^  ■  '^'“[(H,-t)Ps/lw,lf +y 

y  .  _ y 


‘[(H2-t)p;/iw;i]*+y*  ^‘^[(Hj-t)p;/iw;i]Hy" 

Let  Ji.(xi,a)  (i=  1-8)  be  the  asymptotic  parts  of  Ji(xi,a),  as  a  -»«<>,  then  we 


(3.40) 


(3.41) 


have: 


K,3s(xi.t)  =:~r  +  ^2-  Icosatda 


^yj(K+K±hA  r cos(at)e-^'®‘-*‘Ma 

1^14^80  ® 
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Yt(ki%±k)  r  cos(ot)e-”'«-'>do 

’<y»K 

-r.  W,(Hi--Xi)  W3(H,--Xi) 

-x,)f +t"  ™[w,(H,  -X,)]' +t" 


(3.42) 


Using  the  same  technique  we  can  derive: 

t 


^i4s(*i**)  “  Pl5l 


+  Pi6; 


(w,(H,-x,H“+t"  "“[w,(H,-x,)f +t" 

K  (X  tl-o  w^OHj-x,)  w;(H,-x,) 

K23s(x2.t)  - -x,)f  +t^  ^  -x,)]^ + 1* 


(3.43) 

(3.44) 

(3.45) 


K«s(X2,t)  -  ,z  +  -x,)]“  + 

The  coefficients  Xj  (i=85-112)  and  pj  (j=5-20)  are  given  in  Appendix  B. 

The  asymptotic  analysis  for  kernel  KjjXy.t),  K34'(y.t),  K43’(y,t)  and  K44'(y,t)  is  given  in 
Appendix  D. 

3.1.4  Normalization  of  the  singular  integral  equations 

Using  the  following  transformation  formulae: 


Xi=air,  t=aiS,  X2=a2r,  t=a2S. 
y = ^2-bi)r + |(b2+bi),  t = |(b2-bi)s  +  |(b2+b,) 


-l^r^l,  -l^s^l 


then  t+y =|(b2-bi)[s+r+bo), 


-l^r^l,  -l^s^l 


Substituting  above  into  the  four  singular  integral  equations,  we  obtain: 

—  \\  [-^ + Ti&iKi  (r.  s)]<})?  (s)ds  +  f*  ajk^i  (r,  s)(|);  (s)ds 

+  f,  ko  (r.  s)(l);  (s)ds  +  f  k°4  (r.  s)<|);(s)ds  =  -  ^Pt  (r)  (3.46a) 

2  2  •'-1  2Y,4Ey 


£,  (ri  s)<l>i  (s)ds  + Jtajk^  (r,  s)](l>“  (s)ds 
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^i^|;_ki(r.s)«(s)ds+^|;,ki(r,s)$;(s)ds=-^i^^4;«  (3.46b) 


|\,ki(r,s)Jitir(s)ds+£ajki(r,s)K.;(s)ds+|.*;(r)+i»^J^jI4(r,s)(|>;(s)ds 


27f'-»  s-r  s+r+b„  2  •'-»  2E, 

|\,k;,(r,s)](|iJ(s)ds+J^a,ka(r,s)M>;(s)ds+^J^(^^+^|j-^jj^)J5(s)<is 


•^lj^J[,k|.(r,s)(|i;(s)ds-^«(r)+.^i^|‘,k:,(r.s)(t.:(s)ds=^:(r)  (3.46d) 

where  superscript  "o"  is  used  to  denote  the  normalized  quantities. 

For  orthotropic  materials,  the  dominant  parts  of  eqns  (3.46c)  and  (3.46d)  are 
coupled.  In  order  to  solve  the  singular  integral  equations  numerically,  it  is  necessary  to 
decouple  these  dominant  parts  of  eqns  (3.46c)  and  (3.46d).  Using  the  technique  described 
in  [22],  first  let: 


then  eqns  (3.46c)  and  (3.46d)  are  combined  into: 


.  .  B  f  >  (bds 
A<^  +  — 

^  7riJ-»s-r 


+E.  f  ^  i 

jci  J-is+r+bo 


(3.47) 
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Multiplying  both  sides  of  eqn  (3.47)  by  A'*  gives: 

A  +  ±  f‘  A'*B<i>ds  +±  [‘  f‘  A-'Ceds+f  A'‘K<{)ds  =  A“’P  (3.48) 

^  TtiJ-i  s-r  Jti**-*  s+r+bo 


0  -V 


Let  D  =  A-‘B  = 


iJ^  0 


P»  •  d'=a-1b'= 


_ia 

Pi 


then  the  eigenvalues  of  D  can  be  determined  from: 
D-pI  =  0 


giving:  =  =±  =±^ 

Pip4  VPiP4  C 


Now  let  R  be  a  square  matrix  such  that: 
DR  =  RA 


where: 


rr  0 

[o 


(3.49) 


(3.50) 


1  i 


then  we  have  R= 


1 


where  111=— ,  and‘n2=-^ 
Pi  P4 


where 


Introducing :  ^  =  R\|f,  where  y  =  K  ,  then  eqn  (3.48)  is  reduced  to: 

LY  4  J 

1  f.  R-A-B'Ryds^|.^..^-ca^ 

jjjJ-i  s-r  ra*'-!  s+r+bo 
+j‘^R'‘A‘‘K<|)ds  =  R“‘A'’P 

„  _i A  -ip  _  r^i^i  ® 2®2  (*’>  ®)i 

“[aiCj(r,s)  a2C4(r,s)J 


(3.51) 
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‘p,  P4K  p.  P^V^k.! 
P4  Pii^i  P4  P>  - 


(3.52) 


-1  A  -iTrr.  _  t>2  ~k  r^l  (^'»S)  kjCr.S) 


R-‘A'‘KR 


2  Ik'Cr.s)  k 


l(r,s) 

aM. 


_^)+i(  E^+ E.^)  (ML_i!L^)+i  E(^+^) 
>.-b,  V  P4^  Sin.  p/V%P4'  >.  naP4  V%  P.  P4^ 

^.P4  P4^  Vn/p.  P4^  p.  P4  lln.  p.  \n.  P4  J 


(3.53) 


'^■■^■‘=[-::. -tk^rEi'T 

L  \PlP2P4  p4  . 

Finally  eqns  (3.46c)  and  (3.46d)  are  reduced  to: 

Y3(r)+^j!,Vj(s)^-^j!iY5(s)^^j:^^+£iaiCi(r.s)(l)i(s)ds+J_jajC2(r,s^^ 


(3.54) 


b,-b,  f» 


jkj(r,s)vi/|(s)ds+^2Y^£k;(r.s)\K:(s)ds 


1— V  V 

1  yjy^y 

‘  2E. 


.P3(r)+im2— P4(r) 


(3.55) 


hzh-f' 


k; (r, s) yJ (s)ds  +  k“ (r, s)Y4 (s)ds 
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(3.56) 


Note  that  the  functions  Y|(r)  and  yJCr)  are  related  by  the  following  relations: 


V;(r)=iW(r)-i^;(r)] 

V:W=iH^J(r)+«(r)]=-i^^(i) 


(3.57a) 


(3.57b) 


If  both  materials  are  isotropic,  it  is  found  that  =  1,  and  the  above  relations 


reduce  to: 


V:(r)  =  |^i4>S(^)  +  mi  =  -iVsCr) 


(3.57a*) 


(3.46a*) 


where 


(3.57b*) 

Now  to  expressing  eqns  (3.46a, b)  also  in  terms  of  xp^r)  and  have: 

^£jt-^+na,k;(r.s)]((i;(s)ds+f^ajk;’j(r,s)(|i;(s)ds 

+  [T.v;(s)ds  +  T.^(s)lds  =  ip^p?(r) 

2  ^Yl4^y 

f  aik"  (r,s)(}>j  (s)ds  +- f  [— + 7ta2k$i(r,s)](l);(s)ds 
+^f  [T,v;(s)<ls  +  T2^(s)]ds  =  ^;^p;(r) 

T,  =  kJ+i&|i  T,  =  ki-i&l 

T2  =  4+i&i  and  T2  =  ki-i&i 


(3.46b*) 


(3.58abcd) 
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3.1.5  Determination  of  singularities  for  different  crack  geometry 
The  stress  singularity  at  the  crack  tip  varies  with  the  crack  geometry, 
a)  Embedded  cracks  ( |axl<  Hp  |a2l<  H2  and  b2>  0 ) 

For  embedded  cracks,  the  crack  density  functions  are  of  the  following  forms; 


I?(t)  F2(t) 

(i+tra-t)^’  (i+t)“>(i-ty** 


Y3‘(t)  = 


m) 

(i+t)“’(i-ty*‘ ’ 


¥:(t)= 


(i+tra-t/^ 


(3.59abcd) 


where  F°(t)  (i=l,2,3,4)  are  bounded  at  crack  tips,  and  because  of  symmetry  and 


Substituting  eqns  (3.59a,b,c,d)  into  eqns  (3.29a, b,c,d)  and  using  the  following 
formula  [33]: 


.  .  1  f*»<|),(t),  1 

Jf'*i  t-z 


F„(t) 


(t-z,)“-(z2-tr(t-z) 


dt 


F„(zi)cot(7iaJ  F„(z,)cot(7cpJ  ^  ^  _ 

(z2-zi/*(z-z,)“-  (Zj-Zif'CZj-z)***  “ 

where  Gn(z)  is  bounded  at  the  crack  tips. 


Then  eqns  (3.29a,b,c,d)  give: 


F°(-l)cot(7tai)  Fi°(l)cot(7tai)  _  (1-VxyVyx)  .  x 

(2)“'(r+l)“‘  (2r(l-r)“‘  2Yi4Ey 


+  bounded  terms 


(3.60) 


(3.61) 


I^(-l)cot(7ta2)  _  l^(l)cot(7ia2)  ^  (r)  +  bounded  teims  (3.62) 

(2)“‘(r+l)*»  (2)‘^(l-r)“»  2y,4E, 


K(t)  C'^  I^(-l)cot(ita3)  F;(l)cot(7tP3), 

(l+r)“’(l-r)P>  i  (2)^(l+r)“‘  (2)‘^(l-r)^ 

=  miP3(r)  +  im2P4(r)  +  bounded  terms  bj  <  y  <  b2  (3.63) 

fi(r)  C~*  .Fr(-l)cot(7ra4)  F;(l)cot(7rpJ 

(l+r)‘‘^(l-r)'*‘  i'  Util+Tf*  ^ 
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=  -  imiP3(r)  -  m2P4(r)  +  bounded  terms  bi  <  y  <  b2  (3.64) 

From  eqns  (3.61)  and  (3.62),  it  can  be  easily  shown  that: 

(3.65) 

These  are  the  known  results  for  cracks  perpendicular  to  the  interface  [23]. 

Siiuilarly,  from  eqns  (3.63)  and  (3.64)  we  can  derive; 


cot(iia3)-  iC=-i  |  ' 
VP2P3 

COt(7Cp3)  =  iC=  i 

I P2P3 

cot(x04)-  iC-  ij'  '  '  , 
VP2P3 

COt(7Cp4)=  iC-  iJ 

VP2P3 

Using  the  formula: 

cot”‘(z)=-Vlog-^^  (3.66) 

2i  z-i 

We  then  have: 


where 


1  .  Q  1  . 

03  =  - +10),  p3  =  --iCD, 

1  .  o  1  . 

04  =  - -ICO,  P4  =  -+10), 

1  ,  1+^. 

CD  = - log( - -) 

2k  ^  1-C 


(3.67abcd) 


Eqns  (3.67a-d)  are  of  the  same  form  as  those  obtained  for  interface  cracks  between 
two  isotropic  dissimilar  media  [24-29]. 

b)  Matrix  crack  touching  the  interface(a3=H3,b3>0  and  a2<H2) 

In  this  case,  as  Xj  and  t  approach  to  Hj,  the  kernel  kii(xi,t)  becomes  unbounded. 

Let  k,j,(r,s)  be  the  normalized  kernel  of  ki,s(Xi,t),  that  is: 


.0  .  H,(l-s)p5/lw.l+HJwJ 

luU  )  “[H,(l-S)pj/Iw,l+H,lwjl]*-Hf(wir)* 
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.  H,(l-s)Ps/IWtl+H,IWil 

Hi(l-s)Ps/IWil+HilWtl 

"[Hi(l-s)p5/lw,l+H,lw,lf  -Hf  (Wif)* 

.  1  H.(l-s)p5/lw,l+H.lw,l 

«'‘[Hi(l-s)P5/lw3l+Hilw3l]*-Hf(w3r)* 


_  X^IWiI 


^  -+ - - r-5-] 


P5 


^^5  Hi(l-s+(l+r)^)  H,(l-s+(l-r)^) 


Iw  P 
|_r^i2llL' 

Ps 


XgjlWiL  1 _ , _ 1 _ 1 

Ps  Ps 


2P=  ‘H.a-s+(l+r)«)%,(l-s+0-r)!5!fe!) 

Ps  Ps 


\g3lw3l 


2P 


[- 


*  + - - — r-r-l 


‘  H,(l-s+0+r)!^)  H,(l-s+a-r)!^) 

Ps  Ps 


(3,68) 


Substituting  eqn(3.68)  into  eqn(3.29a)  and  moving  the  bounded  tenns  to  the  right  hand 
side  of  the  equation,  we  have: 


-f.'— 

jjj-i  s-r  2pj 


Iw  Iw  1^ 

H,(l-s+(l+r)i-i^)  H,(l-s+(l-r)-^) 
Ps  Ps 


1 


2Ps  H.(l-s+(l+r)'^)  ■  H.(l-s+(l-r)l^) 

Ps  Ps 


-] 
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ti  ^izi^r _ 1[ _ 4. _ _ 

*  9ft  vlw,llw,U^T, . ,,  ^lw,llwjl^ 

^Ps  H,(l-s+(l+r)— Hi(l-s+(l-r)  *  -S 

Ps  Ps 


XT  ^mIWjL  1 


H}<l>r(s)ds 


H,(l-s+(l+r)^)  Hi(l-s+(l-r)^) 

Ps  Ps 

1— V  V 

= - (r)  +  bounded  tenns 

2Y„E, 

Applying  eqns  (3.59a)  and  (3.60)  to  (3.69),  and  setting  ai=Pi=^,  we  have: 


F°(-l)cot(7rY)  _  F.°(l)cot(irY)  X,,  lw,l  F°(l)  1^1. 

(2)'(l+r)»  (2Wl-r)'  2  p, 


,  X,.  Iw.l  F°(l) 


2  Q 

(2)^sin(jry)| 


IWillWgl 


L  P 


s  J 


■[— ^+— 
"(l+r)^  (l-r)’'" 


X„  lw,l  K(X) 


2  D 

*  (2)’ sin  (try  )| 


IwJIWjl 


_[ — 1 1 — ] 
(1+r)^  (l-r)T 


L  Ps 


^  iBl - 

2  P*  ,  w^  ^  (l+r)^  a-r)^ 

(2)^sm(jtY) 

-Ps . 


1  — V  V 

- ^Pi  (r)  +  bounded  terms 


2Yi4Ey 

Multiplying  both  sides  of  eqn  (3.70)  by  (l+r)Y  and  letting  r  -4  -1,  we  have: 


(3.69) 


(3.70) 
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Ps 

‘  -.+x. 


lw,l 

1 

Ps 

IWillWjiT 

Ps  J 

1 

--n 

il 

,ftj 

(3.71) 


Ps  pWillWjlj  Ps 

This  result  is  the  same  as  that  given  in  [23]. 
c)  H-shaped  cracks  (ai=Hi,  bi=0  and  a2>0 ) 

As  aj->  Hi  and  bi-»  0,  the  transverse  crack  intersects  with  the  interface  crack  and  the 
kernels  Ku,  K13,  K14,  K31  and  K41  become  unbounded.  Let  a  be  the  index  of  singularity 
at  the  point  of  intersection,  then  we  have: 

<l‘i(0= — r%ir»  <l>4(s)=-  — 


t“(H,-t)'*^ 


Substituting  above  expressions  into  eqns  (3.29a,c,d),  we  have: 

1  1 


1  f«‘ 


if  ‘  {_J_ 

JtJ-H,  t-Xi 


+  tc^[- 


2P: 


Iw  P  Iw  I' 

S  (H,-t+(H,+x,)-L-)  (H,-t+(H,-x,)-rL-) 

Ps  Ps 


] 


2ft  ‘(H,-.+(H.+x,)!5^)  (H.-t+(H.-x,)™) 

Ps  Ps 


] 


2ft  ‘(H.-t+(H.+x.)M  (H.-t+(H.-x.)^) 

Ps  Ps 


1 _  _ 1 _ X  fi(t) 


2P 


5  (Hi-t+(H,  +  x,) 


Ps 
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+  f^ro  K»(H,-x,)  lw3l(H,-x,)  fiCt) 

Jo  IWil^(Hi-x,)^+t^  lw3l^(Hi-x,)*+t*^“(bi-t)'^ 

It  _ t _ -I  ^4(0  < 

Jo  IWiP  (H,  - Xif  + IWjl^  (H,  - x,)^  + ^  t® (bj  -  tf* 

+  bounded  terms 


Hx-t 


Mu 


H,-t 


l-H,  *  Yi2  !^;^^(H.-t)» 

Ps  Hs 


(H,-t)P,  (H,-I)p. 

Iw,l  lw,l _ P,(t)  .. 

^■P‘(H,-t)^p»  ,  P‘(H.-t)^pj,^^)(H?-iy 


lw,P 


IW3P 


■+y 


[  Pi  A  Cy\4.  SjL  f***  (■—  +  —L.\  ^4(0 

2  2jc  Jo  t+y^  t“  (bj  -  if*  ^ 


1— V  V 

- ^^"^PsCy)  +  bounded  terms 

2E, 


f«‘  _ 

J-H,  ‘2yi37c'-|w,l*y^ 


Iwily 

P5 


Iw3ly 

P5 


+  p. 


■+P10  /TT  _.v2o2  } 


F,(t) 


(H.-tm  3  ■  '^*°(H.-t)^P^  ,^(Hf-tY 

lw,P  ^  IW3P  ^ 


dt 


■I.Pi.  A  (y)+  £2.  f***  (_—  +  -^) - _ dt 

+y  <l»4(y)+  2^  ^  t-y  ^  i+y  t“(b2  - 1)^  ^ 


2G 


■P4(y)  bounded  terms 


»y 


(3.72) 


(3.73) 


(3.74) 


43 


Multiplying  both  sides  of  eqns  (3.72-74)  by  t®  and  letting  t  ->0,  we  then  have  the 
following  three  algebraic  equations  for  Fi(Hi),  F3(0)  and  F4(0)  with  a  yet  to  be 
determined. 


.  cot(7ta)  ^  _ 1 _ ^  _ 1 _ 

(2H,)*  2  pj  (2Hj)“sin(jtaX^r  ^  (2Hi)“an(7ta)[^^^^) 


Xg7  IW3I 


+  _ 1 _ +  _ ! _ 5- 

2  Ps  (2H,)“sin()ta)[i^!^]“  ^  Ps  (2H,)“siii(icaX^I 

Ps  Ps 


K 


■[t^+T^1F3(0)  +  ” 


(7%+7^1F.(0)=0 


2b|cos(— '  IWsf'"’'"  2t^sinS 


(3.75a) 


Pj  lw,l  P5  lw,l  (2H,)“anA 


F,(H,)  +  -£|-F,(0) 


+ 


-^[cot(7ia)+ 


1 

sin(7ta) 


]F4(0)  =  0 


(3.75b) 


-{ 


2Yi3  Iw,I  IwjI  Ps  IwjI 


1^, 

Ps 


(2Hj)“cos(^) 


Since  Fi(Hi),  F3(0)  and  F4(0)  are  generally  assumed  nonzero,  therefore  the 
characteristic  equation  associated  with  eqns  (3.75abc)  has  to  be  zero.  Numerical 
calculations  indicate  that  a  is  0  as  expected. 


44 


3.1.6  Numerical  solutions 

In  this  section,  numerical  procedures  are  outlined  for  solving  the  singular  integral 
equations  for  different  crack  configurations  as  discussed  in  the  previous  section.  The  goal 
is  to  determine  the  stress  intensity  factors  and  the  strain  energy  release  rates  for  each  crack 
configuration,  at  the  crack  tips. 

3.1.6.1  Embedded  cracks 

Rgure  3.1  shows  the  embedded  crack  geometry.  In  this  case,  the  normalized 
singular  integral  equations  take  the  following  form  as  derived  earlier  in  section  3.1.4: 

— f ,  [— + Ttaik'i  (r, s)](l)i  (s)ds  +  (r. s)(l); (s)ds 

+  [*  [XjY3(s)ds  +  Ti^(s)]ds  =  ^-^^^^p?(r)  (3.76a) 

£^aika  (r,s)())"(s)ds ^2J^2(r,s)](l)2(s)ds 

[T2V3(s)ds  +  f2V3(s)]ds  =  — 4^P2(r)  (3.76b) 

ja-'-i  s-r  It  •'-1  s+r+b„  •'-*  •'-* 

+  J’j  kf  (r,  s)\y3  (s)ds + -  k j  (r,  s)v|r;  (s)ds 

2E,  2G^ 

Y4(r) r, V4(s)— - — f , ¥3 (s)  +r. a,C3 (r, (s)ds  +  f*  ajC^ (r,s)(l);(s)ds 

m*'-*  s-r  It  •'-»  s+r+b,  •'-»  •'-1 

+  |‘j  k;  (r.  s)\k;  (s)ds + k;  (r,  s)  v;  (s)ds 


(3.76d) 
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All  the  quantities  are  defined  in  section  3.1.4.  Among  the  four  equations  given 
above,  only  three  are  independent  since  Vj(r)  and  vJCr)  are  related  through  eqn  (3.57b). 

The  first  three  will  be  used. 

Equations  (3.76a)  and  (3.76b)  are  Fredholm  equations  of  first  kind.  Thus  Lobatto- 
Chebydiev  integration  formula  is  used.  Equation  (3.76c)  is  a  Fredholm  equation  of  the 
second  kind,  and  thus  a  Lobatto- Jacobi  quadrature  integration  formula  is  used. 


_  m) 


Vl-s* 


1=0 


here  a  =  -  Pj,  p  =  -  03,  F“(s)  =£  diP/‘^’(s)  and  d,  =  (d^  +  i  d,,) 

1=0 

Using  the  methcxis  described  in  [30,37],  we  have: 

^  w  1 

y  {■^[-i-+7caiK“(ri,Sk)]F“(s,)+  a2WkK2j(ri.s,)F“(Sk)} 
w  Sk-rj 


^  1=0  H 


1  —  V  V 

— ^p?(n) 


(3.77a) 


^  W  1 

£  {a,WkK^(ri,Sk)F“(Sk)  +^[—5—  +  7ta2K^(ri,Sk)]I?(Sk)) 
w  JC  Sk-q 


+ 


bz  ~  bj  ^ 

2  j«o 


£  [T2(ri,Syj)Wyjd,P«‘^\SyP  4*  Ti  (q,  Syj)  W  dj  j^)]  = 
H 


4  *  * 

1-V.yVyx 


P;(n) 


(3.77b) 


Clf  t  C'^f  f  W^d.Pr(s^) 

Jri  ^  s^-r^  It  S  ^  s,j+r^+bo 


+  Z  [aiCi(ri,Sk)If  (Sfc)  +  a2C2(ri,Sk)F2*(Sk)]Wk 
k=l 
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+  •^1  i  rr5(r^.gW„<l.^-^’(s^)+T,V^)W«d,P!'^’(s^^ 


^  1-0  i*l 

(3.77c) 

where  Tj  =  k  *  (ryi,Syj),  T4  =  *^2  (^’yi'Syj) 

and 

(k=l.n) 

(k  =  2,..n-l) 

(k  =  l,...n) 

r7C(2i-l), 

(i  =  l,...n-l) 

2“*»**N![rO+l)rr(a+N+2)  ^  . 

^  r(P+N-i-2)r(a+p+N+3)  ^  ^ 

(yj=0) 

2“^‘N![r(a+l)fr(p+N+2)  ^  . 

r(a+N+2)r(a+p+N+3)  ^  ^ 

(yj  =  N+l) 

2“^*r(a+N+2)r(p+N+2)  1 

”  (N+l)(N+l)!r(a-hp+N-i-3) 

(yj  =  l...N) 

the  abscissa  Syj  are  the  roots  of  the  following  equation; 

a-s,i®)p£;’^"(s«)=o 

and  the  collocation  points  r^j  satisfies  the  equation  of  the  N-th  order  Jacobi  polynomial: 

Wyj  and  p{“^^(Syj)  (1=0, ..N)  are  the  complex  conjugates  of  Wyj  and  P/®^^(Syj)  ieq)ectively. 

At  the  interface,  since  a  and  P  are  complex  numbers,  the  weight  Wyj,  the  abscissa 
Syj  and  the  collocation  point  ryj  will  also  be  complex.  The  technique  of  using  the  complex 
weight,  abscissa  and  the  collocation  points  are  not  new.  One  can  refer  to  [36-39]. 
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Note  that  generally  d(  ( f  =  0,...N)  are  complex  numbers,  and  equation  (3.77c) 
actually  provides  two  sets  of  equations  (for  the  real  and  imaginary  parts).  Together, 
equations  (3.77a-c)  provide  2(n-l)+2(N)  equations.  Since  there  are  2(n)+2(N+l) 
unknowns  be  determined,  the  following  additional  equations  are  required  to  insure  the 
single  valueness  of  the  solution: 

<l)J(s)ds  =  0,  (t>$(s)ds  =  0  and  \|r;(s)ds  =  0 

These  in  turn  provide  four  additional  linear  equations: 

X  Wkii;'(o=o 

k=l 

X  WiF;(g=o 

k=l 


i)  E  W;diP/“^>(s^)  =  0  (3.78abcd) 

1=0  j=l 

Note  also  that  the  last  expression  provides  two  sets  of  equations 

Once  the  values  of  1^(1^),  I^(^)  (k=l,...n)  and  d,  =  dj^  +  id^  ((  =  0,...N)  are 

determined,  the  stress  intensity  factors  can  easily  be  calculated  (see  appendix  E  for  details 
of  the  derivation): 

Defining  the  stress  intensity  factors  as  follow: 
for  transverse  crack  in  matrix  and  fibo” 

at  Xj  “  a^.  k(aj)  “  Trm  ^2(Xj  “  a^)  Ojy(X|,0) 

I|— 

X2  “  ^2*  ~  ^2(X2  ^2^ 

for  interface  cracks: 


aty  =  bi; 


1  1 
P3 


W-i 


VP1P2P3 
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=  ^  (b.  1 


aty  =  b2: 


1  1 
Ps  2G^ 


VP1P2PJ 


=  I^  (b.  -y)“>(y-b.)fc[i^(j,^(H,.y)  -  i  (^i.^^o,„(H„y)  ] 
y-b,  P3  2G,y  VPiPaPs 

After  some  lengthy  calculations  (see  Appendix  E),  we  have 


(3.79ab) 


(3.79c) 


(3  TOU) 

P3  2G^  VPjPzPs 

Jf  the  materials  of  both  layers  are  isotropic,  it  can  be  shown  numerically  that 

1  1  fTT  i-v„v„ 

p3  2G,y  P]P2P3  2Ej 


and  eqns  (3.79c, d)  reduce  to 


k2(bi)-iki(bi)  =2G^p3Vi^F3“(-l)  (3.79c*) 

kjdb)  -  ik,(b,)  =  -2G,P3Vr?I?(l)  (3.79d*) 

The  strain  energy  release  rate  at  the  interface  crack  tip  can  be  derived  as  follows: 

—  =  +  Jl — Ljj^]  (3.80) 

Ay  2E.  p,  ‘  2G,  p,  " 

Again  if  both  layers  are  made  of  isotropic  materials,  the  above  equation  reduces  to 


AE  27t  1  1  „  2  j 


Ay  Vi^2^>Tp3 


(kf  +  kl) 


(3.80*) 
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It  should  be  noted  that  when  upper  and  lower  interface  cracks  meet  (that  is  b|=0), 
eqns  (3.76a-d)  should  be  modified  before  proceeding  to  the  numerical  solution.  Hrst  the 
integral  limits  for  interface  cracks  should  be  extended  fix)m  (0,  to  (-bj,  1)2).  Then,  using 
the  following  relations: 

m = <i>3(-t)  m = •  rM) 

Ki3(xi,t)  =  K53(xi,-t)  Ki,(xi.t)  =  Ki4(xi,-t)  (i  =  1.2) 

Kj3(y,t)  =  Kj3(y,-t)  and  Kj,(y,  t)  =  Kj^Cy.-t)  (j  =  3.4) 
eqns  (3.76a-d)  can  be  reduced  to: 

—  f  [— + jcaikji  (r,  s)](t)°  (s)ds  +  f  (r.  s)<t>2  (s)ds 

jc*'"*  s-r  •'"* 

+  ^|'  [T,v;(s)ds  +  T,^i^(s)]ds  =  i^^^^p?(r)  (3.81a) 

f‘  ajk^  (r,s)<|)i  (s)ds +- f  [— + 7ca2k^(r,s)]<l)2(s)ds 

j-i  jjrf-i  5_r 

+  '|-r,  (T,<rt(s)ds  +  T2ij^(s)lds  =  i^^p;(r)  (3.81b) 

2  2Yi4Ey 

^  J\,c,(r,s)|i;(s)ds+ j\jC,(r,s)iti;(s)ds 
+  -l-f  ,(kt(r,s)X|f;(s) +k;(r,s)v:(s)]ds 


1  — V  V 

laLi)* 


=  ^  ^-paCrl+im;^— P4(r) 

2E,  2G^ 

V“(r)-^f  Y;(s)—  +  f  a,C3(r.s)(|)i  (s)ds+f  a2C4(r,s)(t>$(s)ds 
+  [kj  (r,  s)Y3  (s) + K  (r.  s)  (s)]ds 


(3.81c) 


l-v„v 


=  -nn 


2_L£. 


‘  2E. 


P3(r)-“i^P4(r) 


(3.8  Id) 
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The  above  numerical  procedure  and  the  single  valuedness  conditions  can  also  be 
applied  to  solve  eqns  (3.81a-c). 

3.1.6.2  Matrix  cracks  touching  the  interface  (a^sHp  bi>0) 

When  matrix  cracks  touch  the  interface,  the  kernel  I^j(x,,t)  becomes  unbounded. 
And  the  singular  behavior  at  the  crack  tip  changes.  Now,  redefining  (j)® (s)  as: 


K(s) 

a-sy 


the  power  of  singularity  y  can  be  determined  from  eqn  (3.7 1). 

Since  y  is  generally  different  than  0.5,  in  the  numerical  procedures  described  for 
eqns  (3.78a-c),  we  now  use  the  Lobatto-Jacobi  quadrature  for  integrals  containing  <l>“(s). 
The  rest  of  the  definitions  and  numerical  calculations  in  eqns  (3.78a-c)  remain  the  same. 
The  modified  version  takes  the  following  form; 


t,  {^[-ri-7+ra,K“s(r;,s;))ir(sLHW;a,K;^r:.s;)ir(5,)+ajW,K*(^,s,)IJ(,i,)) 

K  Sk-rj 


1=0  j=l 


^Yw^y 


a  W  1 

£  {a,W;K5,(ri,s;)F“(s;)-»-^[-i-  +xa2K5,(ri,Sk)]F“^ 

-  Sk-q 


K 


(3.82a) 


+ 


by-b, 

2 


I’dO 


X  [T2(ri,Syj)Wyj4P;‘^\Syj)  +  Tj  (rj.Syj)  W,j  dl  (Syj)] 

H 


1  *  * 
2y;4e; 


vlir,) 


(3.82b) 


y  w^d^pr-P^(s^)  A  w,id.pr^\s^) 

TCI  ls=0  Jssl  Sjj  ““  Tyi  TC  1^:0  jj-l  Sjj  +  Tjj  +  bjj 


k*l 


[aiCi(ryi,sJI?^(Sk)W‘  +a2C2(ryi,s^ 
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+  i  [T,(r^.gW„4fJ“'>(s,p  +  T,(r^,s^W«d,p!“'’(s,p] 

^  l«0  j»l 


=  mi^-^^P3(r,i)+im2^p,(r^) 


(3.82c) 


where 


r(n+2-2Y)  ^ 


Wk 


,  2^'^~^(n-l)![r(l-Y)f _ l_ 


(k=l,n) 

(k  =  2,..n-l) 


r(n+2-2Y)  [PiT'^^Cs;)]* 

the  abscissa  s'^  are  the  roots  of  the  following  equation: 

and  the  collocation  pcmts  rl  are  the  roots  of  the  following  equation 
Pili‘-^'‘'(r;)  =  0 

The  methods  of  evaluating  the  weight,  abscissa  Wj^,  t,j  for  Lobatto-Chebyshev 
integration  rule,  Wj,  Syj  for  Lobatto-Jacobi  integration  rule  and  their  collocation  points  q 
and  ryj  will  be  the  same  as  those  described  in  section  3. 1 .6. 1 . 

The  single  valuedness  conditions  become: 

k«l 

2  WtIJ(y  =  0 

k^l 

N  n 


I E  Wjd,pr>(s,j)=o 

1=0  H 


(3.83abcd) 


The  expressions  for  stress  intensity  factors  and  strain  energy  release  rates  given  in 
the  previous  section  can  still  be  used  in  this  case  except  for  k(ai) 

Redefining  k(aj)  =  Lim  2^(x2+H2)ra2y(x2,0) 


we  have  (see  Appendix  E): 
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k(ai) 


Kl/Ps  ,  .  k|/p5 

(l-v;v;)sin(Try)^^-Bw:H/p,f  ‘“O-^hl/PsF 


kl/Ps 

[Kjwil/pjf 


+x, 
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kl/Ps 

[kkl/Psf 


] 


(3.84) 


3.1.6.3  H-shaped  cracks  (ajsHj,  bi=0) 

When  Ej  =  Hi  and  bj  =  0,  the  transverse  matrix  cracks  intersect  with  the  interface 
crack  to  form  H-shaped  cracks  as  shown  in  Fig.  3.2.  This  configuration  has  been  observed 
during  the  tensile  test  of  the  uniaxial  ceramic  matrix  composites.  Since  the  power  of 
singularity  at  the  point  of  intersection  is  zero,  which  has  been  confirmed  numerically,  (j)®  (s) 
will  no  longer  be  singular  at  this  crack  tip.  Also,  the  kernels  Kn,  Kjj,  K14,  K31  and  K41 
become  unbounded.  In  this  case,  because  bj  =  0,  the  numerical  calculations  will  be  based 
on  eqns  (3.81a-c)  instead  of  eqns  (3.76a-c). 

Hrst,  the  singular  parts  of  the  kernels  Kjj,  K^,  Kj^,  Kjj  and  K41  in  eqns  (3.81a-c) 
will  be  extracted  and  integrated  in  closed  form.  Then  the  nwdified  equations  take  the 
following  form: 

-f,  [— +7taik“s(r,s)](|)®(s)ds+ ai[k“  (r,s)-k“  ,(r.s)](l)r(s)ds 

+  £_  ajkJ(r,s)$;(s)ds  +  -^fi  (T„v;(s)ds  +  T„N^(s)]ds 

+  -^Jl,  [fr,-T, J v;(s)ds  +  (T|  -  T„)^(s)]ds  =  (3.85a) 

aik^(r,s)(l)i(s)ds  +  ;^£j  [•^+7ta2k^(r,s)](l);s)ds 

+^J‘  rr,v;(s)ds+T.;K(s)]ds=^^^p;(r) 

/  zy,4iiy 


(3.85b) 
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V3(r)+^  f  VJCs)—  +  f,  a,c,,(r,s)(l)f (s)ds  +  f‘  ai[Ci(r,s)  -  c,,(r,s)](l>i(s)ds 

jjj  v-i  g_j  v-i  rf-i 

+  £_  ajCj(r,sH;(s)ds  +  -|-||jk;(r,s)\|/;(s)+k;(r,s)v;(s)]ds 

=  °°i  '■”-Pi(0+iPiT;^P4(r)  (3-85C) 

2E^  2G^ 

where  T,.  =  kS.  +  i^k;.  Tu=ki.->^K. 

and  c,.(r,s)  =  [iii  +  i^Jll 

Pi  Pi  Vni 

Now  defining  <|)f  (s)  =  I^(s),  the  definitions  for  the  rest  of  the  displacement  density 
functions  remain  the  same  as  in  section  3.I.6.I.  Then  we  have: 


{^[ 

Jt 


Sk-q 


+  7caiK«  .(r;,s;)]F“(s;)  +  aiW^K"  (ri.s,)  -  K“  ,(r:.s;)]F“(s;) 


ajWtKi(r|,Si)IJ(Sk))+-^|;  X  tT„(ri.s,j)W„4^“®(s^)  +  T,.(ri.s^)Wnd,p!“®(s^)] 


1=0  j=I 


b  ^  “ 


X  X  fT,(ri,s^T,.(r|.s,j))W„<i,I^<“«(s,j)  +  [T,{r„s,j)-T„(ri.s^]WndiPl“'’(s^l 


1=0  H 


2YmE, 


(3.77a) 


X  (a,w;Ki(r:,s;)F,'-(4)  ■3-[— 5—  +  ra^Kifr^SiMF^CSi)) 

It 


k=l 

N  n 


+  i  [T2(r,SypWyjd,P.^“«(SyP  +  T2(r,s^)W,id.Pi“^\syj)]  =  i^ 

^  1=0  H  -^llA^y 


A  W^diP'-^’Cs^)  C"f  f  WndiP^’CSa) 
Jci^  ^  s^-r^  JC  S  ^  s,j+r^+b^ 


(3.77b) 
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+  S  W'k  a,Ci,(r;  .sJF,- (4)  +  £  w;  ai[  c,(r:  .<)  -  c„(r:  ,4)1^ (4) 


k«l 


k«l 


+ 1,  Wia,Cj(rA)^(Sk)At  t  [T,(r^^W„4,Ff"'>(s,p 

k«i  ^  1=0  H 


+  T,(r^,S^W«d,Pi““(s,pi  =  m,^--:^p,(r^)+imj;^p,(r^)  (3.77c) 


l-v^v, 


2E 


2G, 


where 


w;  = 
w:  = 


n(n-l) 

2 


1 


(k=l.n) 

(k=2,..n-l) 


n(n-l)[P„_,(4)f 

the  abscissa  are  the  roots  of  the  following  equation: 

(l-4')Pn^’(s‘k)  =  0 

and  s’j.  are  the  weights  and  abscissa  of  the  Lobatto-Legendre  integration  rule. 
The  collocation  points  q  are  the  roots  of  the  following  equation 

pir^(r;)=o 


Again,  the  methods  of  evaluating  the  weight  and  abscissa  Wj^,  t|^  for  Lobatto- 
Chebyshev  integration  rule;  Wj,  Syj  for  Lobatto-Jacobi  integration  rule  and  then- 
collocation  points  q  and  ryj  remain  the  same  as  those  given  in  section  3. 1.6.1. 

The  remaining  four  additional  equations  take  the  following  form: 


Fx“(-1)  =  0 

(3.86a) 

f;(1)=o 

(3.86b) 

<|)2(s)ds  =  0 

(3.86c) 

and 

|^<|)i(s)ds  =  ||^(|)3(s)ds 

(3.86d) 

The  first  two  equations  originate  fiom  the  requirement  that  <|)“(s)  remains  finite  at 
both  crack  dps.  The  third  equation  is  deduced  from  the  single  valuedness  condition  for 
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f  cracks  in  the  fiber  if  82  is  not  zero.  And  the  last  equation  states  the  condition  of 

continuity  of  displacement  at  the  point  of  intersection. 

Eqns  (3.86ab)  can  be  directly  used  in  the  numerical  calculations.  Eqn  (3.86c)  can 
be  written  as: 


Z  Wkl^(tk)  =  0  (3.86c*) 

Using  the  fact  that 

J^(i)3(s)ds=-|-J‘^ 

=  •^•7]!,  V3(s)ds 

=  .i^do - - - , 

2  ^cosh(7tco) 


we  have 


J_^(}),(s)ds  =-i^do 


K 

cosh(7tco) 


Since  I  <l)i(s)ds  is  real,  therefore, 

J-Hi 

Re(do)  =  0  (3.86d*) 

This  result  gives  the  last  equation  that  is  needed  for  the  numerical  calculations. 

The  calculation  of  the  stress  intensity  factors  at  the  embedded  crack  tips  remain  the 
same,  and  needs  no  further  elaboration. 

3.1.7  Results  and  discussion  for  singular  integration  formulation 

It  is  well  known  that  the  interface  bonding  strength  contributes  to  the  toughness  of 
ceramic  matrix  composites.  A  "weak"  bonding  strength  generally  improves  the  toughness 
of  the  composite.  A  direct  indication  of  weak  bonding  between  the  matrix  and  the  fiber 
interface  is  the  existence  of  interface  crack  or  defects  either  due  to  manufacturing  or  due 
to  debonding  during  loading.  In  the  experiment,  it  was  observed  that  under  tensile  loading. 
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matrix  cracks  generated  first  and  then  propagated  to  the  interface.  Once  the  matrix  cracks 
reached  the  interface,  they  did  not  penetrate  through  the  fibers.  This  means  the  firacture 
energy  available  at  the  crack  tip  is  less  than  the  firacture  energy  of  the  fiber,  and  is  larger 
than  the  firacture  energy  along  the  interface  due  the  weak  interface.  Thus,  the  energy  is 
released  along  the  interface  causing  interface  debonding  or  generating  interface  cracks. 
The  actual  role  that  the  interface  plays  has  drawn  a  lot  of  attention  by  both  involved  in 
research  and  manufacturing  of  ceramic  composites.  One  can  refer  to  [45,46]  for  detailed 
discussions.  The  goal  here  is  to  explain  or  simulate  the  damage  mechanism  of  ceramic 
matrix  composites  by  analyzing  the  interactions  between  matrix  and  interface  cracks  using 
different  crack  configurations. 

The  numerical  calculation  for  embedded  crack  configuration  has  been  carried  out 
first.  To  check  the  correcmess  of  the  numerical  procedure,  some  simple  crack 
configurations  have  been  studied  first  For  periodic  transverse  cracks  in  the  matrix  and  the 
fibers  without  interface  cracks,  the  results  match  those  given  in  [23,44];  For  interface 
cracks  without  transverse  cracks,  if  both  Hj  and  are  set  very  large  and  the  materials  of 
both  layers  are  isotropic,  the  problem  is  reduced  to  an  interface  crack  in  an  infinite  domain 
and  the  results  fix)m  [25,43]  for  stress  intensity  factor  at  the  interface  were  recovered. 
Table  3.1  gives  the  results  of  stress  intensity  factors  at  the  interface  for  four  different 
materials  combinations.  As  can  be  seen  fi'om  the  table,  the  two  results  are  almost  identicaL 
Also,  the  strain  energy  release  rates  at  the  interface  for  the  above  four  material 
combinations  were  calculated  using  eqn  (3.80*).  For  isotropic  materials,  the  exact 
expression  of  the  strain  energy  release  rate  at  the  interface  was  first  given  in  [40]  as 
follows: 

AE  TC  (1^2  ,  1^2) 

Ay  2  ^ll^l2[(l+K,)^l2  +  (l+K2)^ll]  ‘  ^ 

where  |ij  and  iXj  are  the  shear  moduli  of  the  fiber  and  the  matrix  respectively. 


(3.87) 
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Table  3.1  Comparison  of  stress  intensity  factors 


Materials 

K.(b) 

K.(b) 

CoVb 

CToVb 

Layer  1 

Layer  2 

Results  from 

f251 

eqn  (3.79d) 

[251 

eqn  (3.79d) 

Aluminum 

ISlISflilll 

1 

1.0000 

-0.1342 

-0.13421 

Steel 

1 

1.0000 

-0.1443 

-0.14431 

Steel 

Aluminum 

1 

1.0000 

-0.09158 

-0.091576 

bficalon 

CASH 

1 

1.0000 

-0.071076 

-0.071074 

The  results  using  eqn  (3.80*)  differed  with  that  given  in  eqn  (3.87)  by  a  factor  of 
VJ .  The  strain  energy  release  rates  using  both  eqn  (3.80*)  and  eqn  (3.87)  are  compared  in 
Table  3.2.  The  error  is  within  3  percent.  Considering  the  fact  that  all  the  derivations  are 
based  on  the  orthotropic  materials,  and  the  isotropic  behavior  was  approximated  in  the 
calculation,  the  results  are  quite  good.  Table  3.3  gives  the  Young's  moduli  and  the 
poisson's  ratios  for  the  materials  used  in  the  calculation. 


Table  3.2  Comparison  of  strain  energy  release  rates 


Materials 

Strain  energy  release  rates 

Layer  1 

Layer  2 

2  AE  .  O-TN 
- m  (3.87) 

7t  Ay 

V2-— in(3.80*) 

K  Ay 

Aluminum 

0.397614x10-* 

0.402399x10-* 

Steel 

0.384227x10-* 

0.389954x10-* 

Steel 

Aluminum 

0.239708x10-^ 

0.239753x10-* 

Nicalon 

CASH 

0.203262x10-* 

0.208256x10-* 
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Table  3.3  Material  constants 


Material 

Young's  modulus 

Poisson's  ratio 

E  Ob  in. -2) 

Aluminum 

1x10’ 

0.3 

Steel 

3x10’ 

0.3 

Epoxy 

4.5x10* 

0.35 

Nicalon 

2.83x10’ 

0.28 

CASH 

13.2x10’ 

0.28 

Next,  the  effects  of  the  transverse  and  interface  cracks  on  the  tensile  damage 
behavior  of  the  ceramic  matrix  composite  are  investigated.  Suppose  there  exist  both  small 
transverse  matrix  cracks  and  interface  cracks  or  defects.  First,  let  the  matrix  cracks  be 
allowed  to  propagate  and  the  interface  crack  be  fixed  at  a  small  value  b/H2=0.01.  The 
stress  intensi^  factors  at  both  the  matrix  crack  and  the  interface  crack  tips  are  calculated 
and  plotted  against  the  increasing  matrix  crack  length.  Hg.  3.3  shows  the  mode  I  stress 
intensity  factors  at  the  matrix  crack  tips  versus  the  matrix  crack  length  for  both  Vf  =  30% 
and  Vf  =  40%.  The  stress  intensity  factors  increased  first  and  then  decreased  as  the  matrix 
crack  approached  the  interface.  The  decrease  of  the  stress  intensity  factors  as  the  matrix 
crack  approached  the  interface  is  because  of  the  fibCT  constraint  effect  The  mode  I  stress 
intensity  factors  at  the  interface,  however,  remained  unchanged  at  first  and  increased 
rapidly  as  shown  in  Fig.  3.4  when  the  matrix  crack  approached  the  interface.  Fig.  3.5  gives 
the  strain  energy  release  rate  at  the  interface  versus  the  matrix  crack  length  for  two 
different  fiber  volume  fractions.  The  influence  of  fiber  volume  fractions  can  also  be  noted 
here.  Next,  let  the  matrix  crack  length  be  fixed  (&2=0-91i^  and  the  interface  crack  length 
be  allowed  to  increase,  then  the  stress  intensity  factors  at  the  matrix  crack  tips  and 
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interface  crack  tips  were  calculated  and  plotted  in  Fig.  3.6  and  Fig.  3.7.  respectively.  The 
mode  I  stress  intensity  factors  at  the  transverse  matrix  crack  tips  exhibited  a  steady 
increase  until  the  interface  crack  length  reached  a  certain  value  (b=H2).  Fig.  3.8  shows  the 
strain  energy  release  rate  at  the  interface  versus  the  interface  crack  length.  The  following 
conclusion  can  be  drawn  fiom  the  above  two  cases:  If  both  matrix  crack  and  interface 
crack  exist  in  the  composite,  the  matrix  crack  will  first  propagate  because  of  the  increase 
of  stress  intensity  factors  at  the  matrix  crack  tips  and  the  low  fracture  toughness  of  the 
matrix.  Then  the  crack  will  be  arrested  when  it  reaches  the  interface  because  of  the  fiber 
constraint  effect.  On  the  other  hand,  if  the  fracture  toughness  along  the  interface  is  small, 
the  interface  crack  will  grow  due  to  an  increase  in  the  strain  energy  release  rate  at  the 
interface.  Once  the  interface  crack  or  debonding  length  reaches  a  certain  value,  the  fibw 
constraint  effect  is  almost  lost  and  the  stress  intensity  factors  at  the  matrix  crack  tips  start 
to  grow  again  as  can  be  seen  from  Fig.  3.9  where  b/H2  =  0.55.  Eventually  the  matrix  crack 
will  reach  the  interface  boundary.  This  has  been  observed  during  testing.  Another 
important  role  that  the  interface  crack  plays,  is  to  prevent  the  matrix  cracks  from 
extending  self-similarly  into  the  fibas.  As  can  be  seen  fiom  the  test  results,  fibers  do  not 
break  at  the  same  locations  as  the  matrix  resulting  in  longer  fiber  pull-out  length.  The 
variations  of  the  stress  intensity  factors  versus  the  increase  of  matrix  crack  length  for  two 
different  values  of  interface  crack  length  are  shown  in  Rgure  3.9a  and  Figure  3.9b  for 
fiber  volume  fractions  of  30%  and  40%  respectively. 
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3.2  Finite  Element  analysis 

The  H-shaped  crack  configuration  is  adopted  and  used  to  study  the  effects  of 
interface  cracking  and  its  progression  on  the  tensile  behavior  of  ceramic  matrix 
con^sites. 

3.2.1  Description  of  the  H-shaped  crack  model 

From  the  experiments,  it  is  observed  that,  after  the  appearance  of  initial  transverse 
matrix  cracks,  the  stress-strain  behavior  starts  to  become  nonlinear.  >^th  certain  increase 
of  the  tensile  load  (about  10  ksi),  a  regularly  spaced  multiple  transverse  matrix  crack 
pattern  is  formed.  Once  this  stage  is  reached,  further  increase  of  load  generates  few  new 
matrix  cracks  until  final  failure.  Therefore,  during  the  period  fiom  which  saturated  matrix 
crack  formed  to  eventual  failure,  interface  debonding,  matrix  crack  opening  and  sequential 
fiber  breaking  dominate  the  failure  process.  To  capture  the  effect  of  damage  accumulation 
due  to  interface  debonding,  a  configuration  of  periodic  H-cracks  as  shown  in  Hg.  3.2  is 
used  where  the  interface  crack  length  is  allowed  to  increase  with  increasing  load.  The 
popular  ABAQUS  finite  element  code  is  used  with  strain  energy  release  rate  as  the 
criterion  for  the  interface  crack  propagation.  Fust,  the  strain  energy  was  calculated  for  a 
certain  interface  crack  length,  then  the  interface  crack  length  is  allowed  to  increase  by 
releasing  the  node  at  the  crack  tip  and  the  strain  energy  is  again  calculated.  Thus  an 
estimation  of  the  strain  energy  release  rate  G  at  the  interface  crack  tip  can  be  determined 
using  the  following  expression: 

G  =  ^n-yn±i  (n=0,1.2, . ) 

tAD 

where  Uq  and  Un+i  are  the  strain  energies  associated  with  the  two  subsequent  crack 
lengths,  AD  is  the  interface  crack  increment  and  t  is  the  thickness  of  the  specimen.  If  this 
G  value  is  greater  than  G^  (critical  strain  energy  release  rate),  then  the  crack  is  assumed 
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to  propagate.  In  practice,  the  Gcr  value  at  the  interface  is  extremely  difficult  to  measure. 
Here  we  estimate  0^  by  matching  one  point  of  the  stress-strain  curve  with  the  calculated 
value.  Table  3.4  shows  the  computed  G^r  for  both  Vf  =  30%  and  Vf  =  40%. 


Table  3.4  Estimated  critical  strain  energy  release  rate 


Fiber  volume  fraction 

L/2H, 

Estimated  Gcr 

30% 

2.5 

40% 

3.0 

62.70  Obs/in) 

3.2.2  Finite  element  results  and  discussion 

Before  the  finite  element  calculation  is  carried  out,  one  has  to  determine  the 
matrix  crack  spacing  L  and  the  fiber  spacing  2H2Used  in  the  model  as  shown  in  Fig.  3.10. 
The  ratio  affects  the  results  significantly.  Matrix  crack  spacing  is  determined  by 
averaging  the  matrix  crack  length  fiom  the  micrographs  taken  during  the  erqjeriment  It  is 
found  that  the  matrix  crack  spacing  decreases  with  increasing  fiber  volume  firaction  of  the 
specimen.  The  average  matrix  crack  spacing  for  Vf  =  30%  was  found  to  be  L  =  166  pni 
And  for  Vf  =  40%,  the  value  is  L  =  111  pm.  The  fiber  spacing  varies  with  fiber  volume 
firaction,  the  fiber  diameter  and  how  evenly  the  fiber  filament  is  distributed  in  the  matrix 
media.  The  average  fiber  diameter  for  Nicalon  fibCT  is  about  16  pm  as  measured  from  the 
test  results.  Consider  the  case  where  fibers  are  evenly  distributed  in  the  matrix  media, 
then  for  Vf  =  30%,  2H2  =  (0.7)(16)/0.3  =  37.3  pm  and  for  Vf  =  40%,  2H2=  (0.6)(16)/0.4 
=  24  pm.  Therefore  we  have: 

—  =  1 66/37  =  4.44  for  Vf  =  30% 

2PI2 

and  •;—  =  1 1 1/24  =  4.625  for  Vf  =  40% 


62 


Considering  the  fact  that  fibers  are  randomly  distributed  in  the  matrix  and  that  the 
fiber  diameter  is  not  constant,  the  above  ratios  were  varied  between  2.0  to  3.0  in  the 
calculations.  Hg  3.1 1  is  a  comparison  of  finite  clement  results  with  the  test  results  for  Vf 
=  30%.  Fig.  3.12  shows  the  comparison  for  Vf  =  40  %.  It  can  be  seen  that  finite  element 
results  match  the  test  results  quite  well  at  the  start  of  the  non-Unear  stress-strain  curve. 
Then  the  two  results  deviate  somewhat  as  the  load  is  furthw  increased.  This  is  because  at 
higher  loading  the  fibers  start  to  break  as  can  be  seen  from  the  micrographs  taken  during 
the  tests.  This  contributes  to  the  flattering  of  the  actual  stress-strain  curve.  It  may  also  be 
noted  that  the  finite  element  results  for  Vf  =  30%  match  the  stress-strain  diagram  better 
than  those  for  Vf  =  40%.  This  can  be  explained  as  follows:  with  more  fibers  in  the  40% 
fiber  volume  fiiaction  specimen,  more  sequential  fiber  breakings  occurred  during  the  final 
stage  of  the  experiment  after  the  saturated  multiple  matrix  crack  pattern  was  formed.  This 
can  not  be  simulated  by  considering  interface  cracks  only.  Fig.  3.13  gives  the  stress 
distribudon  in  both  fiber  and  matrix  near  the  interface  crack  tips.  Fig.  3.14  gives  the 
normal  stress  distribudon  along  the  interface 
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Fig.  2.1  Experimental  setup 


Fig.  2.2  Specimen  geometry  ( Unit  in  inches) 


Fig.  2.4  Tensile  tests  on  Nicalon/CAS  n 
with  Vj=30%  at  room  temperature 


to 


(IS5I)  SS9J}g 


Strain 


Matrix  crack  initiated  Fully  developed  matrix  crack 

away  from  the  defects 

Fig.  2.6  Another  pattern  of  matrix  crack  initiation. 


Fig.  2.7  Micrographs  showing  the  location  of  the  first  matrix 
cracks  initiation  and  their  development. 


(C)  at  55  ksi 


(D)  at  20  ksi 


Fig.  2.8  Tensile  damage  of  Nicalon/CAS  n  specimen  with  Vf=40% 
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Load  =  60  lbs 


Load  =  49.6  lbs  Load  =  33  lbs 

Fig.  2.9  Further  loading  of  the  specimen  is  possible 

after  failure  for  thin  NIcalon/CAS  n  specimen. 


Location  and  the  damage  pattern  just  before 
major  fracture  developed. 


Major  fracture  on  the  specimen 


Fig.  2.10  Brittle  behavior  of  the  thick  Nicalon/CAS  n  specimens. 


Fig.  2.11  Long  fiber  full-out  length  for  Nicalon/CAS  H  with  Vf  =  40% 


Fig.  2.12  Short  fiber  puli-out  iength  for  Nicalon/CAS  n  with  Vf=30%. 
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Fig.  2.13  Very  smooth  surfaces  on  the  fibers  being  pull>out 


Fig.  2.15  Typical  damage  pattern  for  NIcalon/CAS  n  specimen 
taken  at  different  magnifications:  (a)  300x;  (b)  800x;  (c)  3000x. 
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Fig.  2.16  Typical  fiber  breaking  pattern  for  Nicalon/CAS  n  specimen, 


Fig.  2.17  Radom  fiber  breaking  on  Nicalon/CAS  n  specimen 


(c)  36.5  ksi  (d)  51.5  ksi 

Fig.  2.18  Matrix  cracks  density  vs.  the  tensile  stress 
for  Nicalon/CAS  n  specimen  with  Vf=30% 


Fig.  2.19  Matrix  cracks  density  vs.  the  tensile  stress  at  room 
temperature  for  Nicalon/CAS  n  specimen  with  Vf=40% 


Matrix  crack  density  vs.  stress  for  Nicalon/ CAS  II 
specimens  with  Vf=30%  at  room  temperature 


Stress  (ksi) 


Room  temperature 


T  =  250  OC 


T  =  400  OC  T  =  600 


Fig.  2.26  Final  failure  surfaces  at  different  temperatures 
for  Nicalon/CAS  n  with  Vf  =  30%  . 


Fig.  2.27  Development  of  matrix  cracks  at  T=600®C 
for  Nicalon/CAS  II  specimen  with  Vf=30%. 


(a)  Pre-existing  crack  (b)  21.1  ksi 

before  loading 


(c)  32.5  ksi  (d)  40.8  ksi 

Fig.  2.29  Development  of  matrix  cracks  at  T=400  OC 
for  Nicalon/CAS  n  specimen  with  Vf=40% 


(c)  56.1  ksi 


(d)  64.1  ksi 


Fig.  2.31  Progression  of  fiber  crack  opening  at  T=400  OC 
for  Nicalon/CAS  n  specimen  with  Vf=40% 


(b)  56.5  ksi  (c)  64.1  ksi 

Fig.  2.32  Progression  of  fiber  crack  opening  at  T=600 
for  Nicalon/CAS  n  specimen  with  Vf=30% 


Fig.  2.33  Tensile  Stren^  of  Nicalon/CASn 

with  V^0%  at  various  temperatures 


Temperature  (®C) 


g.  3.3  Stress  intensity  factors  at  the  matrix  crack 

for  Nicalon/ CAS  II  (b/H2=0.001) 


Matrix  crack  length  (a2/H2) 


Fig.  3.5  Strain  energy  release  rate  at  the  interface 
for  Nicalon/CAS  II  with  (b/H2=0.001) 


Matrix  crack  length  (aa/Ha) 


3.6  Stress  intensity  factors  at  the  matrix 

for  Nicalon/CAS  II  (a2/H2=0.9) 


Interface  crack  length  (b/H2) 


.7  Stress  intensity  factors  at  the  interface 

for  Nicalon/ CAS  II  (a2/H2=0.9) 


Interface  crack  length  (b/H2) 


Matrix  crack  length  {o.z/'Rz) 
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Vo 


Fig.  3.10  Finite  element  model 


Fig.  3.11  Tensile  tests  on  Nicalon/CAS  H 

with  Vj^30%  at  room  temperature 
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(IS51)  SS3JJS 


Strain 


(XS5[)  SS94S 


.1^  Normal  stress  distribution  at  the  interface 
th  different  interf0Ci0l  crack  length  CL/2H^=2.5D 


[Tsd)  ssaj>s 


element  number  along  the  interface 


Appendix  A 

From  eqns(1.9ab),  one  obtains; 


_  _  i.  r  [A(a)sinh(WiaHj)  +  C(a)sinh(w3aHj)]asmayda 

»y  jc  Jo 

—  f  f  +  '^''(^^)e"*'''3^®y/P5]sinatsinotHjdtda  (A.1) 

^Jo  Jo  Yi3  Ps  Ps  112 

^^2(~^2»y)  _  £  f  [A*(a)sinh(wi*(xH2)  +  C(a)sinh(w3*0(H2)]asmayda 
oy  Jt  Jo 

2  r  p  il!2(^.i^eMwi*lay/P5*+i|3^i^-lw3*lay/p5*] 

JcJo  Jo  Yi3  Ps  Ps  Yi2 

^i^i»y)^  —  f  [p7A(a)cosh(wi(xHi)  +  pgC(a)cosh(w3(xH,)]acosayda 
It  Jo 

-fir 

+ ^?^ioSign(w3)e-lw3lay/p5]sinatcos(xHidtda  (A.3) 

Ps  Yi2 

^2(-H2.yL  1  f  [^*A*(a)cosh(w,*aH2)  +  p8*C*(a)cosh(w3*aH2)]otcosayda 
^y  It  jq 

-  irr 

+  *  "p  j0*sign(w3*)e’'w3**®^y/P5*]sinatoos(xH2dtda  (A.4) 

P5*Yi2* 

Using  the  following  integral  formulas: 

J  e'^*sinbtdt=^^  and  j*  e'^^osbt  dt  =  ^2+^  ^>0 


(A.5ab) 


equations  (A.1-A.4)  reduce  to: 


au.(H„y) 

dy 


7t  Jo 


[A(a)sinh(w,(xHi)  +  C(a)sinh(w3(xHi)]asinayda 

^»2(-^>yLl  r  [A*(a)sinh(w,*aH2)  +  C*(a)sinh(w3*ocH2)]asinayda 
^y  JC  Jo 


ii 


where 


— i.  f  [B,A(a)cosh(wiaHi)  +  PgC(a)cosh(w3aHi)]cccosayda 
oy  It  Jo 

-ff 

^^('^•yL  £  r  [p^*A*(a)cosh(wi*(xH2)  +  pg*C*(a)cosh(w3*aH2)](xcosayda 
oy  7C  Jo 

2  r» 


■If  $^Ay)d. 


(A.4*) 


I  Wily 

IWlir _ &5 _ 

iw,ly 

IWgl  Ps 

~p5  (iw:^!y)!. 


jwjk 


My 

A 


■1 


■] 


!wi*!y. 
_§s!_ 


lwi*ly 

p5* 


Iwa*ly  iw»*ly 

lw,*l,  Ps*  _  P.*  , 

- 2“^* - Slit 


■] 


i7 


-] 


P,nSign(w,)lw,hrii,  H.-t  H.+t 

-  fe  ;;;w,(H,.t)r^.cH.V 


.  Po*sign(w,*)lw,*l 
54(t,y)  =  '^’  ^ 


i(wi*)lw,*l  H,-t  _  H,+t 


P,n*sign(w,*)lw,*hr,  1  *  H,-t 


Ps* 


^(lw,*ly)^ 


P5*2 


(A.6a) 


(A.6b) 


(A.6c) 


(A.6d) 


Applying  (3.14ab)  and  (3.21cd)  to  eqns(A.l*)-(A.4*); 

4>4(y)  =  -  —  f  [A(a)sinh(WiaHj)  +  C(a)sinh(w3aHi) 

Jt  Jo 

+  A*(a)sinh(Wi*aH2)  +  C*(a)sinh(w3*(xH2)](xsinayda 

r 

(|)j(y)  =  —  f  [^A(a)cosh(wiaHi)  +  P8C(a)cosh(w3(xHi) 

7t  Jq 

-  p7*A*(a)cosh(w3*aH2)  -  pg*C*(a)cosh(w3*aH2)]acosayda 

Taking  inverse  cosine  and  sine  Fourier  transforms  for  both  sides  of  eqn(A.7a)  and 
(A.7b)  respectively,  one  gets: 


A(a)sinh(WiaHi)  +  C(a)sinh(w3CtHi) 

+  A*(a)sinh(w,*aH2)  +  C*(a)sinh(w3*(xH2)  =  Ri(a)  (A.8a) 

and 

p7A(a)cosh(wiCxHj)  +  p8C(a)cosh(w3aHi) 

-  P7*A*(a)cosh(wi*cxH2)  -  pg*C*(a)cosh(w3*aH2)  =  R2(oc)  (A.8b) 

where 

1’4(y)staaydy  +  ;|J^  ^.(t.y)staayduly 
•^(“>=^1  -IS(y)co»y<Iy+^||  ^,(t.y)cosaydtdy 

r  <A.9b) 


Using  the  following  integrals: 


iv 


J  Gi(t,y)sinaydy  =  - 

+ —  e'®(^i'*'*)P5/^'^3b  (A.  10a) 

J  G2(t,y)smaydy  =  - ^e*®^032‘OP5*/'wi*l  -  e-a(H2+t)P5*/lwj*I) 

+ ?^e-a(H2-t)p5*/lw3*l  -  e-a(H2+t)p5*/lw3*l)  (A.iOb) 

2Yi2* 

J*  G3(t,y)cosaydy  =  ^ign(Wj)P9(c‘®^OHrOP5/'wil  -  e-aOH|+t)P5/IW|l) 

-  ^ign(w3)P  jo^\e‘®(HrOP5/'w3l  _  e-a(Hi+t)P5/lw3l)  (A.  10c) 

^  112 

J  G4(t,y)cosaydy  =  ^ign(Wj*)P9*(e‘®0H2‘*)p5*/*W|*l  _  e-a(H2+t)P5*/IW|*l 

-  |sign(w3*)pio*^e-a(H2-t)p5*/Jw3*l  -  e-a(H2+t)p5*yiw3*l)  (a.  lOd) 

eqns(A.9a)  and  (A.9b)  are  finally  reduced  to: 


wb. 


<l)4(y)sinaydy 


^4(XY  I  +^^e■®^CHl-t)P5/lw3l](Jt 

‘V-ai 


Yi2 


f  <|i2(t)[-e’“(^'*)p5*/lwi*l  +  ^e-a(H2-t)P5*/lw3*l  jdt 


and 


(A.lla) 


R2(a) 


=-Lf 

2a  I 
Jh, 


<t>3(y)cosaydy 


<|)i(t)[sign(Wi)p9e-«(Hi-0p5/'Wil  -  sign(W3)pi(^-«(Hi-t)p5/lw3l  ]dt 


I  ‘{*2(0[sign(Wi*)p9*e-C‘(H2-t)p5*/Iwj*l 


-  Sign(W3*)Pio''^^^'®^2'^^P5*^W3*^  ]dt 
Tl2 


(A.llb) 


Applying  the  continuity  conditions  (3.1  Sab)  and  (3.14ab),  one  derives  the 
following: 


I®*  2[y3A(a)cosh(WiaHj)  +  Y4C(a)cosh(w3aHj) 

Jo 

-XjY3*A*(a)cosh(Wj*otH2)  -  X,Y4*C*(a)cosh(w3*otH2)]occosayda 


=s-  r®®  f^i^^sinatdt[Yie"'^i'®y/P5  -  Y^^'^^j'^y/PslcosocHida 

JoJo 

I-  f“  f“i  ^‘^^^atdt[Y,*e-l’*'i*l“y/P5*  -  Y,»^-IW3*l“y/p5*]cosaHjda 

JoJo 


(A.12a) 


and 


I  ®®  2[Y9A(a)sinh(WiaHi)  +  YioC(a)sinh(w3aH,) 

Jo 

+  X2Y9*A*(a)sinh(W|*(xH2)  +  X2Yio’‘'C*(a)sinh(w3*aH2)]asinayda 
=-  r®®  r^i^^sinatdtYij[e‘*Wi'“y/Ps  -  e’l'^s'^y/PsisinaHida 

JoJo 

-  r~  f»i  ^^4^atdtY„*[e-l»'i*l“y/Pi*  -  e-lws'^y/Ps'lsinoHida  (A.12W 

JoJo 

Applying  inverve  Fourier  cosine  and  sine  transforms  to  (A.12a,b)  respectively, 


we  have: 

[Y3A(a)cosh(w3aHj)  +  Y4C(a)cosh(w3aHi) 


(A.13a) 


vi 


-  X,,Y3*A*(a)cosh(W|*aH2)  -  XiY4*C*(a)cosh(w3*(xH2)]a= RjCa) 
and 

[Y9A(a)sinh(WiOcHi)  +  YioC(a)sinh(w3aH|) 

+  X2Y9*A*(a)sinh(W|*aH2)  +  X2Yio*C*(a)sinh(w3*(xH2)]a= R4(a) 
where 


(A.13b) 


R3(a)=-“ 
2 


iff 


^.fcy)co.ayd<dy+fj~J>i  ^ 


^^^^^^20*y)cosaydtdy  (A.  14a) 


R4(a)=^"J~J®i^^3Msinaydtdy-^J‘”J^i  ^^^^4(t,y)sinaydtdy  (A.14b) 


2  /•„  To 
'0-^0 


and 


- Hiit 


P5^ 


+(Hrt)2 


— -  .  Y  - iilil - .1 


•] 


+  [Yi 


S  —  f—  V  * - ^  Y  \ - ^2.1^ — 1 


+  [Yi*- 


P5« 

H?+t - _  ^  Yufl^ - Kj+t - 


■] 


Iwily  Iw^ly 

lw,ly  lw,ly 


■] 


lw,*ly 

P.* 


Iw,*ly 

id. 


•] 


-Yii*[: 


lwt*ly 

-iiL. 


Ps*^ 


lw3*!y 

P.* 


(Iw,*ly)^^ 
P5* 


(A.15a) 


(A.  15b) 


(A.  15c) 


(A.15d) 


vii 


With  the  help  of  the  following  integrals: 


f-S.(t.y)cosaydy  -  - -  e-«(H.«)|5*dw,l) 

Jo  ‘ 

EMii§s^g.a(Hi-t)p5/lw3l  _  (A.  16a) 

2Yi2Iw/ 

f~S,(..,)cosaydy  = 

Jo  ‘ 

^7.*Yi  1  -  e'C^2+0P5’'‘/*'^3’''b  (A.  16b) 

2  712*1^3**  ' 

f~S,(t,y)stoaydy  =  -  e-a(H,«)P^»,l) 

Jo  ^ 

_  2E]fuPs/e-a(Hj-t)p5/lw3l  _  e-aCHi+OPs/lWsI)  (A.  16c) 

2 IW3I  ^ 

K(W)  W  =  .fu^.-a(H.-Ofe*/lw,*l  -  e-a(H,«)fe*/'w,*l) 

Jo  ' 

_Hu:§s*  a(H,.t)ps*/lwj»l  _  e-o(Hj+t)p5»/IWj»l)  (A.16d) 

2  IW3*I 

R3(a)  and  R4(a)  are  finally  reduced  to: 


(A.  17a) 


J  •JL* 


(A.17b) 
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Appendix  B 

Definition  of  material  constants  for  both  layers.  A  superscript  ♦  will  be  used  for  the 
material  in  second  layer. 

Yj=1+2V^  1  +  Y3  =  W,+Vxyp7  Y4  =  W3+VyxPg 

Ps  Ps 

Y5  =  Vxy+^  Y6=Vxy  +  ^^  Y7  =  VxyWi  +  p7  Y8  =  Vxy  W3  +  pg 

|w| 

Y9 = -1  +  P7W,  Yio  =  -1  +  Ms  '^11  =  “  R  +sign(Wi)p<, 

Ps 


Yi2  =  -y  +sign(w3)Pio 

Ps 

Yi3  =  sign(wi)p9-^ 

Ti2 

sign(w3)Pio 

V  -iri-Is-IiLi 

Y»yJ 

1  e;  i-v„v„ 

‘-E,  i-v;v; 

^  G, 

GO. 

It 

^4“^lY3p7“P7Y3 

X5 = Xi  y4P7~P$Y3 

^  =  Y9-Yio 

X7=Y9-Y9X2 

Xg  rs  Y9  “  Y10X2 

^10  X^ 

^11  “^Y9 

^12“  ^6^4 

Xi3  = 

^14  =  ^eYs 

^15  “  ^6^7 

^16  =  ^^7 

^17  “^P?  ■*'^4p8 

^18  ==  ^4^ 

^19  +^)P8 

^0  =  ^5P7 

^  "^3p8 

^='Y3p7 

P7P8 

X24  =  ^10^16  ~^^16 

^25  =  ^10^17 

X26  =  X12XJ5  -  XjoXjg 

^27  =  X9X20  ~  ^13^16 

^g  =  ^^19 

X29  =  X12X19  -  X13X17 

X30  =  Xj3Xjg  -  X12X20 

X31  = 

X32  =  XiiXj7 

^3  =  ^16^12  “  ^iAi8 

X34  =  X9X22  -  ^14^16 

X35  =  X9X21 

Xi6  —  ~^14^17  ~  ^Ai2 

X37  =  X,i4Xig  -  X22Xi2 

Xgg  =  XjjXjg  -  P7P7X9 

Xy9=y^% 

^40=  ^15^17  ~  ^12^ 

^41  “  P7P7^12  “  ^15^18 

^42~^^16 

X43  =  X3X17 

^44  “  ^^18 

IX 


^45  •"  ^10^16  ~  ^11^16 
^48  “  ^14^16  *"  ^10^ 
x^i = ^14X19 + y^x^i 

X54  =  -XjoX23 

X57  =  X3Xj9 

^60  “  79(^4^19  “  ^5^17) 
^3~“^(73P|  +^5) 


^46  ~  ^11^0  “  ^13^16 
^49  “  ^13^  "*  ^14^0 
^2  “  ^IoPjPt  “  ^15^16 

^55  ”  ^13^  "■  ^15^19 
=  “  ^4^16(^8  “  Y9) 
^61  “73^7(^10“^) 
X|54  =  -  P7p7(A.io  -  X9) 


X47  -  — XnXi9 

^53  =  ^15^0  “ 

X55  =  -  A3X20 

X59  =  —  X5Xj5(Y9  —  X7) 

^2  “  ^sCTj^H  +  ^4^1) 

j  =  Xg(A4A23  “  ^7^17) 


^<6  ~  ^(P7^19  “  ^5^)  ^7  ~  “  ^4^16 

X59  =  A5A17  -  X4AJ9 

A71  =  ^(P»  (Kll  ~  ^)  p8^(79  *"  ^)) 

^3  “  ^1(^10  ~ 

A75  =  A9A5  —  ^{3X3  +  Pg  (A9Y3  —  Aj4A3) 

X77  =  P7  (A10P7  -  A15A3) 

A79  =  X3(A4Pi  —  A5P7  ) 


A^g  —  ^5^16 

A70  ~  ^(P?  ^10  ^11)  Pg^4^8  “  79)) 

A72  ”  p7  (^11^  ■"  ^^13)  ■*■  Pi  (^^12  ~  ^11^4) 
A74  =  A,j2A3  —  ^10X4  +  P7  (X14A3  —  AjoYs) 

X76  =  X23(Aio  —  A9) 

X7g  =  P«  (X15A3  -  P7A9) 

Xgo  =  X24  +  A25  +  A26  +  A27  -  Ajg  +  A29 


Xgi = -  X70  -A7i  -A72  +  sign(wi)p9(A73  +  A74  +  A75) 


+  "f^CXTe  +  ^7  +  ■'['^^(A-43  -  A57  +  X79) 

H  hi 

A82  ~  ^58  —  ^9  “  ^60  +  Sign(Wj)P9(Agi  +  A42  +  Ag3) 

+  ■?^(X54  +  X65  +  Agg  )  -  (  Ag7  +  Agg  +  Ag9) 

|w,|  hi 


Ag3  =  •^(A7o  +  A71  +  A72)  -  J^sign(w3)p,o(  A73  +  A74  +  A75) 

7i2  7u 

_  TilPjTi^  ^  ^  ^  ^  j  ^  T”^|*  (A43  —  A57  +  A79) 
7i2  IWjI  IWjI 

Aj4  =  "^(Agg  +  Aj9  +  Ago)  —  ■^sign(w3)Pio(Agi  +  Ag2 + Ag3) 

7i2  7i2 


_  ^  +  Xg5  +  X^)  +  Y (X-CT  +  +  ^69) 

Y«  lw,|  |wj| 


^-Y7^1  \  ■  ^  =  Y8^2„-^-j:"  ^7“Y7^3  j,  ^8-Y7^  j, 

Yi3Ym'''I0  TaTu^M  TaTw^M  TaTu^ 

^89  =  -  X70  -  -  X72  -  sign(w’)pj  (^^3  +  ^74  +  X75) 

■“  +  ^7  +  ^8)  “  ^Yn  “  ^57  +  ^9) 

^0 = -  ^58  -  ^9  -  ^60  -  sign(w[)p  J  (X^i  +  X^2 + ^63) 

“  "**  +  ^66)  ~^Yll  ^  ^ 

*  • 

X91  =  -^(X^O  +  +  ^2  )  +  ^Sign(W3)P,o  (X73  +  ^4  +  X75) 

Yi2  Ya 

+  ^iY23*'T%(^6  +  Vz  +  ^s)  +  ^2Yu  T%(^43  “  ^  +  ^9) 

Yi2  KI  [WjI 

•  * 

X92  =  "^^(^58  ^59  ^6o)  ■*■  ■^^sign(Wj)PjQ  (Xgj  +  X42  +  ^63) 

Ya  Ya 

^=Y7^9  .  5,  ■  ^  =Y8^-7-  Y“  ^5=Y7^1-r— Y“  ^=Y8^~ 

TaTi4'''80  ian4'''80  lalM^so 

X97  =  -X43  —  ^4^  —  X47  +  sign(Wj)P9(X48  +  X.50  +  X,5i) 


YaYi4^8 


+  ■^^■(X>52  +  X54  +  X55)  -  ^”^1*  (X,42  +  X55  +  X57) 

nil  nil 

^8  “  "  ^1  “  ^2  “  ^3  ■*■  Sign(Wi)P9(X34  +  X35  +  X35) 

^  +  ^9  +  ^4o)  ■  ■’■  ^43  +  ^44) 


XI 


^  “  ^^<5  +  ^4«  +  ^47)  “  ^48  ^50  +  ^5l) 

Yu  Yu 

•  ■*■  ^54  ■*■  ^55)  ■*■  (^42  +  ^  +  ^) 

Yu  ral  IwjI 

Xjoo  =  •^(^1  +  ^2  ■*■  ^3)  “  ~^sign(w3)Pio(  X34  +  X35  +  Xgg) 

Yu  Yu 


■?^^8  +  ^9  +  ^4o)  +  l^*%(“^42"'^43  +  ^44) 

Yu  IW3I 


^01  -  ^Y7„  ^102  -  ^8Yt„ 

YuYm^w  YuYm^m 

^103  =  ^9Y7  „  ^104  =  ^100 Y«  ""”•"7 

YuYm^w  YuYm^so 

Xj05  =  ~  X,45  —  ^4^  —  X47  —  sign( Wj  )P  9  ^4g  +  X50  +  X5|) 

“  Yl  ]^(^52  +  ^54  +  ^55)  ~^Yu  j^^42  +  ^  +  ^57) 

^106  “  ~  ^1  ”"  ^2  “■  ^3  ”  sign( Wj  )P9  (X^  +  X35  +  Xgg) 


-  +  ^9  +  ^4o)  ~^Ya  “  ^43  +  ^44) 

*  » 

^107  “  "^^(^45  ^46  ■**  ^47)  ■*■  ■^^sign(Wj)P,Q  ^4g  +  X5Q  +  X,5j) 

Yu  Yu 

+  ^1  Y2‘^T^(^52  +  ^54  +  ^55)  ■*■  ^Yu  T^(^42  +  ^56  +  ^57) 
Yu  IwjI  {WjI 

•  • 

^108  “  ■*■  ^2  ^3)  +  ■^sign(w3)PjQ  (X34  +  X^s  +  Xgg) 

Yu  Yu 

+  ^  Y2-% j^(^8  +  ^9  +  ^4o)  +  ^zYu  “  ^43  +  ^44) 


Appendix  C 


f,(a)  =  tanh(w*,aH2)tanh(w*3aH2)  fjCa)  =  tanh(w*3aH2)tanh(wiaHi) 
f3(a)  =  tanh(w*3aH2)tanh(w3aH,)  f4(a)  =  tanh(w*,aH2)tanh(w3aHi) 
fjCa)  =  tanh(w*i(XH2)taiih(Wi<xHi)  f^(a)  -  tanh(wjaHi)tanh(w3aH,) 
f((X)  =  X24fi(tt)  +  ^9^6^®) 

gi(a)  =  X7otanh(w*3aH2)  +  X7itanh(w*i(XH2)  +  X72tanh(w3aHi) 
h|((X)  =  ^3f|((x)  +  X74f3(cc)  +  A75f4((X) 
nijCct)  =  A^gfiCct)  +  Xjjf^(fx)  +  X7gf4(ct) 
ni(a)  =  X43tanh(w*3aH2)  -  X37tanh(w*iOcH2)  +  X79tanh(w3CXHi) 
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Appendix  D 

Consider  the  following  integrals  in  eqns(1.29c,d): 


I 


,1)2  /•bj 

K33’(y,t)4>3(t)dt 


|.b. 

Kj,'(y,t)Wt)dt 

Jb, 


I*!’!  rh 

K43'(y.t)(l)3(t)dt  K44’(y.t)<l)4(t)dt 

h^  Jb, 


(D.labcd) 


where 


1  r  ch(w,(xx,)  chCwottXi) 

K„’(y.t)  =  um  - 

Jo 

..  if  .ch(Wiaxj)Y,h,(a)  ch(w,(xx,)  Y4h,(a), 

=  lim  ”  I  [“77 — „  .  V  .  +  .  /  ^  TT » ''t7^  x"]cosaycosatda 

..-.H,  jt  I  ch(WjOtH,)  f(a)  ch(w3aHj)  f(a)  ^ 

Jo 

U(w,aHi)  f(a)  chCwjOH,)  f(a)  y)+«’««>+yM<>“ 

Jo 

w=  to  if 

Jo 

-..ir 

•'O 

_  -L  rrCh(wjaxj) Y3gj(a)  chfwattx,)^ 

— WiOH.)  f(a)  -"chCwsOH.)  f(a)  ][™a(t-y)+sma(t+y)]da 

Jo 

. .  if  ,sh(w,ax,).  ,  ,  sh(w,(xx,) . 

K„(y.O  =  hm  -  ^  )J,3(y.a)  +  ch(w 

Jo 


,ch(w,ax,),  ,  ,  ch(w,ax,) . 


cli(w,ax,)Y,s,(a) 

^chCwiCtH,)  f(a)  chCwjOcHi)  f(a)  J«>saysmataa 


(D.2a) 


(D.2b) 


xix 


sh(w,ax,)Yoh,(a) 
U(w,cxHi)  f(a) 


+ 


sh(w.,ax.)Yinh,(a) 

chCwsOHi)  f(a) 


]sinaycosatda 


sh(w,axi)Yoh,(a) 
xh(wiaHj)  f(a) 


sh(w,ax,)Yinh,(a) 
chCwjCxHi)  f(a) 


][-sma(t-y)+sina(t+y)]da 


(D.2c) 


K,4-(y,t)  =  lim 


i 


sh(w,ax,) 

^ch(wiaHj) 


Ji5(y.“)  + 


sh(w,axt) 

ch(w3aHi) 


JjgCy.a)  jsinatda 


1  r  ,sh(w,ax,)Ygg,(a)  shCw^ax,)  YiQg2(«l 

%  I  ch(wiaHj)  f(a)  chCwjOcH,)  f(a)  ^ 

Jo 

(D2.d) 

i^27c  I  ^ch(WiaHi)  f(a)  chCwaOHi)  f(a)  ^ 

Jo 


As  a  approaches  infinty,  we  have: 


=  ^ie-«W|(Hi-x,)  +  ^  e-aw3(Hi-Xi)  (D.Sa) 

'^chCwiOHi)  f(a)  ch(w3aHi)  f(a) 

ch(w,c«,)K,(a)  ^  ch(w,o«,)m(a)  ^  ^,e-«w,(H.-x.)  (ojb) 

xh(wiCxHi)  f(a)  ch(w3CxHi)  f(a) 

A(w,ax,)^  ^  |,e-aw.CH,-x.)  +  (d.3c) 

‘^chCwjOHi)  f(a)  ch(w3aHi)  f(a)  ^ 

^  e-awi(Hi-Xi)  +  (D.3d) 

xh(W|0cHi)  f(a)  ch(w3aHi)  f(a) 


where  ^re  material  constants  and  are  expressed  as: 

1  ^24  +  +  ^26  +X27  ■"  X2g  +X29 

K  Y4(Xfii+Xfi2+Xfi,) _ 

^  X24  +  X25  +  X25  +X27  “  Xjg  ■*'X29 

t  _ _ Tt(X7n'*'X7l+X7?) _ 

3  X24  +  X25  +  x^6  +X27  -  X28  +X29 

E  _ _ YaCXss+XijQ+X^n) _ 

X24  X25  +  X26  +X27  “  X28  +X29 

t  _ _ Y9(X73+X74+X>7g) _ 

^  X24  +  X25  +  X25  +X27  -  X28  +X29 


(D.4a) 

(D.4b) 

(D.4c) 

(D.4d) 

(D.4e) 


XX 


^4  ^  ^  ~  ^  +^ 

- Y9(^0'^^I+^?) _ 

^4  ■*■  ^  ^  +^7  ~  ^ 

_ Ylo(^?i!t+^<io+^fin) _ 

^4  ■*■  ^  ^  +^7  “  ^  +^9 


Using  the  following  formulas: 


(D.4f) 

(D4.g) 

(D4.h) 


and 


to  I  e‘®''^i(Hi-Xi)[cosa(t-y)  +  cosa(t+y)3da = Jt[5(t-y)  +  5(t+y)] 

Jo 

^  j*  ■•■  sina(t+y)]da  =  (~ + 


(D5.a) 


(D.5b) 


The  asymptotic  values  of  eqns(D.2abcd)  become: 


^  2%  I  +  ^2®'^^3(HrXi)][cosa(t-y)  +  cosa(t+y)]da 

Jo 


=  ■^[5(t-y)  +  8(t+y)] 


(D.6a) 


""  27C  I  +  44®'®^'^30Hi-Xi)][sina(t-y)  +  sina(t+y)]da 

Jo 


27cVy^t+y^ 


(D.6b) 


^  27t  I  +  ^e'®W3(H,-Xi)][_sin(X(t-y)  +  sina(t+y)]da 

Jo 


=  £l(_-L+X) 
27c^  t-y^tV 


(D.6c) 


-  to 


2jj  J*  [^7®"®"^!^*'^!^  +  ^ge‘®^'''3(Hi-Xi)][cosa(t-y)  -  cosa(t+y)]da 


xxu 


J  K,4'(y.«»4(t)<it=-^Wy)  + J  K„(y.i»,(t)dt 


where 


KsjCy.t) 


••>=“? if  ft 

Jo 


ch(w,ax,)Y^h,(a)  ch(w,ax,)  Y.h,(a) 
chCwiOH,)  f(a)  chCwjOH,)  f(a)  ^ 


•  Kie‘“'''i(HrXi)  +  ^2e“®^w3(HrXi)]}[cosa(t-y)  +  cosa(t+y)]da 


cxx,)Y,g,(a)  ch(w,cxx,)Y^g,(a) 
)  f(a)  chCwjOHi)  f(a)  ^ 


„  ,  V  ..  if  ,,ch(w,cxx,, 

K34(y,t)=  -  lim  r”  I  { 

^  *,-.11.  lit  I  ‘  chCwiOHj 

•^0 

•  +  ^4e"®^30Hr*i)]}[sina(t-y)  +  sma(t+y)]da 

rr  15m  f  ,rSh(w,ax,)Yght(a)  .  sh(w,ax,)Yinhi(a) 

™27t  ^U(wi(xHi)  f(a)  ■^chCwaOH.)  f(a)  ^ 

Jo 

■  +  ^e‘^W3(^i"^i)]}[-sina(t-y)  +  sina(t+y)]da 

K  fvrt=  lim-i-T  .  .sKwaOx,)  Ying2(«) 

^271 1  ^^ch(w,ccHi)  f(a)  cMwoOH.)  f(a) 

-  K7e‘“WiCHi-Xi)  +  ^ge’®^w3(Hi-Xi)])[cosa(t-y)  -  cosa(t+y)]da 


(D.8d) 


(D.9a) 


(D.9b) 


(D.9c) 


(D.9d) 
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Appendix  E  Derivation  of  stress  intensity  factors  and  strain  energy  release  rates 
i)  Case  of  embedded  cracks 


The  stress  intensity  factors  at  the  transverse  crack  tips  are  defined  as: 
atxi  =  a,; 


k(ai)  =  IimV2(Xi-ai)aiy(Xi.O) 


at  X2  =  82: 


k(a2)  =  Lim 


From  eqns  (3,29a),  for  Ixil  >  aj  we  obtain: 

where  <Jr,(x„0)  is  a  bounded  function,  and 

F,(l)  _ 

^  ,^ai-V  (t-ai)‘"(t+ai)‘'* 
ftom  [33]  if  we  introduce  the  sectionally  holomoiphic  function  (p(z)  as 


dien  we  have: 


^  Jl*'-'  t-z 


Fi(ai) 


+  <Po(z) 


(2a,y"(z+a,)‘^  '  (2a,)‘"(z-ai)‘'" 

Substituting  (E.3)  into  (E,2)  and  from  the  definition  of  (E.la),  we  have: 


Similaily 


k(a.)  =  -  iF,(a,)  =  - 

X(l~V^Vyj)  ^^2 


Now  consider  the  stress  intensity  factors  on  the  interface. 
From  Eqn  (3.55),  we  obtain: 


I 


P1P2P3 


p3  2G^ 


(E.la) 

(E.lb) 

(E.2) 


(E.3) 


lrl>l 


xxiv 


where  = 

Again  from  [33],  we  obtain: 


1  j:  .  F,(r)  ] 

where  Fjj(r)  is  the  principal  part  of  vJ(r)  at  infinity  and  is  bounded  TTierefore: 


-I?(r)+F.(r)] 


l-v„v„ 


V  P1P2P3 


P^'2,G^ 


(r+l)“*(r-iy 


Multiply  both  sides  of  the  above  equation  by  -i,  then  we  have: 


Now  if  we  define  the  following  stress  intensity  factors  at  the  interface: 

— -^kjCr)  -  i  (r) 

VP1P2P3  2E, 

=  .^(r.l)”.(r-l)^l±^a„(H.,)  -  ] 

Aen  we  obtain: 

In  order  to  derive  Ae  etqtression  for  Ae  strain  energy  release  rate  at  Ae  mteiface 
crack  tip,  we  first  express  Y|(r)  m  terms  of  k,  and  kj ,  Aen  from  (E.4),  we  have 

■■7Rr'ii*'-'E^‘'<STr3^7  "<■  <"> 

Noting  also  Aat: 


.l-v„v„ 


xxvi 


X 


[cos(CDlog(^))  +  i  sin((olog(^))] 


(E.8*) 


The  oscillatory  behavior  in  (E.8*)  can  be  removed  using  the  techniques  described 
in  [40,42].  Considering  the  fact  that  the  stress  intensity  factors  and  the  displacement 
derivatives  just  away  firom  the  interface  crack  tips  should  have  definite  values.  TTierefore 
we  can  chose  T|  and  q  in  such  a  way, 


cologC^)  =  0  and  ci)log(^)  s  0 

It  It 

In  fact,  for  most  material  combinations,  that  above  two  quantities  are  indeed  very  small. 
Then  from  eqns  (E.7)  and  (E.8*),  we  have 

Av  =(v,-V2)  s  ^k2 

VI-C2  P3 


and  Au  =  (Uj-Uj)  =  a 

Assume  the  crack  firont  was  closed  along  the  interface  by  an  amount  Aa 

then  using  the  method  given  in  [44],  we  obtain  the  increment  of  the  strain  energy: 

ff+At 

AE  =  2j  [a„Au(r-Aa)  +  <y^v(r-Aa)]dr 


= ^  r  1  1  1 

p2  2E,  Pj  2G^ 

2E. 

Therefore,  the  strain  energy  release  rate  will  be: 

AH  2g  ,  1  l-v„v„  ■ ;  1  1  ,  2, 

Aa  P2  2E.  ‘  P2  2G,  ^' 


2  j*r+Ai  Vl+  Aa—r 

Vr-1 

k^]|Aa 


dr 


(E.9) 


Fbr  Isotropic  materials  we  have: 


xxvu 


^  =  L.(  icJ+kSl 

ii)  Case  of  transverse  crack  touching  the  interface 
Define  the  stress  intensity  factor 

k(ai)=  Lim  2^(x2+H2)yo^(x2,0) 

from  eqn  (3.29b),  we  have 

°  ■  f  K„,(x,.t)(l.,(t)dt  +  <^,(X2,0) 

where  a$y(x2,0)  is  bounded  function  and 


2jc  1  1 


9i(t) 


F.(t)  F.(t)e^ 


Again  define: 
then  we  have: 


r^dt 

JtJ-h,  t-z 


<p(z)= 


F.(-H,)e^ 


F.(H,) 


+  <Po(z) 


(E.9*) 


(E.IO) 


(E.11) 


(2Hi)Tsin7ry(z+H,)’f  (2Hi)Tsin7ry(z-Hi)’ 

Substituting  (E,12)  back  into  (E.  1 1),  and  then  from  the  definition  of  (E.IO),  we  have: 

I  ,  (H,)^r;E;F,°q)  kl/P,  '  kl/fe 

*  (l-v>^)sin(7ry)  ‘°‘|]w;|w,|/P5]’  ‘*”Dw3*|w,|/ps]^ 


(E.12) 


+  X. 


103 


K|/ps 

[KKl/psf 


+  X- 


101 


K|/p5  3 

[|wj(w;|/P5y 


(E.13) 


The  stress  intensity  factors  at  other  crack  tips  and  the  strain  energy  release  rate  at 
the  interface  crack  tip  remain  the  same. 


iii)  H-shaped  cracks 

The  derivation  of  stress  intensity  factors  and  of  the  strain  energy  release  rates  at 
the  interface  crack  tip  remain  the  same  as  those  given  in  section  i). 


